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Abstract 

We show that the Green's function of a two dimensional fcrmion with a 
modified dispersion relation and short distance parameter a is given by the 
Lerch zeta function. The Green's function is defined on a cylinder of radius 
R and we show that the condition R — a yields the Riemann zeta function as 
a quantum transition amplitude for the fermion. We formulate the Riemann 
hypothesis physically as a nonzero condition on the transition amplitude be- 
tween two special states associated with the point of origin and a point half way 
around the cylinder each of which are fixed points of a Z2 transformation. By 
studying partial sums we show that that the transition amplitude formulation is 
analogous to neutrino mixing in a low dimensional context. We also derive the 
thermal partition function of the fermionic theory and the thermal divergence 
at temperature I /a. In an alternative harmonic oscillator formalism we discuss 
the relation to the fermionic description of two dimensional string theory and 
matrix models. Finally we derive various representations of the Green's func- 
tion using energy momentum integrals, point particle path integrals, and string 
propagators. 



1 Introduction 

The Riemann Hypothesis (that the Riemann zeta function has nontrivial zeros only 
when the real part of its argument is one half) embodies one of the deepest mysteries 
in mathematics. Mathematics and physics have similar origins and many seemingly 
abstract results in mathematics turn out to have physical realizations. Because of the 
depth of the Riemann hypothesis one may hope not only that a physical realization 
of the hypothesis may lead to its solution, but also that its formulation in physical 
terms may uncover physical ideas perhaps related to short distance physics and 
string theory. There are several examples of a physical realization of mathematical 
developments leading to physical advances. The relation of the Atiyah-Singer index 
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theorem to anomalies, instantons, and chiral fermions is one example, the relation 
of Chern-Simons gauge theory and knot theory to three manifold theory is another, 
also one has the Frenkel-Kac vertex operator construction of infinite Lie algebras 
and the relation to the heterotic string. Thus one can expect that the physical 
realization of mathematical problems like the Riemann hypothesis has benefit for 
both mathematics and physics. 

The search for a physical realization of the Riemann Hypothesis has a long his- 
tory going back to Hilbert and Polya who looked for a quantum hermitian operator 
with eigenvalues at the Riemann zeros. More recent approaches include the relation 
of the Riemann zeros to a Jost function in scattering theory [1], the relation of the 
zeta function to quantum chaos [2], and the relation of the Riemann hypothesis to 
noncommutative geometry [3]. More information about recent approaches to the 
Riemann hypothesis are described in [4]. 

In this paper this we give a new physical realization of the Riemann hypothesis 

in terms of a 1 + 1 dimensional fermionic theory whose short distance dynamics are 

modified in a definite way. The short distance dynamics introduce a new parameter 

of dimension length or inverse energy. Within this theory we derive the propagator or 

Green's function which gives the transition amplitude to observe the fermion at any 

point in space-time (x, t) given that its initial position and time. When the spatial 

coordinate x livetakes its values on a circle S 1 of circumference L = 2irR we show 

that this Green's function is proportional to the Lerch zeta function with arguments 
. x it R. ^ , , , L R . 

( , — , — ) . Only when the ratio = — = 1 , that is when the radius of 

2t:R a a zira a 

the universe equals the new short distance scale , do we obtain the Riemann zeta 

function, so in this sense the Riemann zeta function probes short distance behavior of 

the Green's function. By considering initial states of the fermionic particle that are 

statistically distributed, (essentially Boltzmann distributed with respect to energy 

and a mixing parameter a ) we show that the transition amplitude is modified by the 

x %t (j R 

Lerch zeta function with parameters ( , — I — , — ) . The Riemann hypothesis is 

2irR a a a 

L 

formulated physically as saying that the transition probability to a position — half 

way across the 1 + 1 dimensional universe is strictly nonzero for all times and values 

of the mixing parameter in the range — < — < 1 . 

2 a 

This paper is organized as follows. In section 2 we give a review of certain Dirich- 
let series and zeta functions that will occur in this paper. In section 3 we discuss 
the transition amplitudes and Green's function of the 1+1 dimensional fermionic 
theory in the first quantized and second quantized field theory point of view and 
explain the relation to the zeta function. In section 4 we introduce the notion of 
a Boltzmann distributed initial state and introduce the mixing parameter a . We 
formulate the Riemann hypothesis as saying that the transition rate is never zero 
for the transition from initial state parameterized by a and x = to a final position 
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state at x = — = ttR . We discuss a simplified analysis in terms of partial sums and 
show how the behavior of the transition probabilities is reflected in the phenomena 
of neutrino mixing. We also discuss the a parameter from several points of view 
including the integral over theta functions, target space duality, and space-time re- 
flections. In section 5 we discuss the statistical mechanics of the nonstandard 1+1 
dimensional fermionic theory and its relation to multiplicative number theory and 
thermal divergence. In section 6 we discuss an alternate formulation in terms a non- 
standard oscillator and discuss the relation to the fermionic form of two dimensional 
string theory and matrix models. In section 7 we discuss the relationship between 
ground ring operators, algebraic curves, modular forms and L-series. In section 8 
we discuss other representations of the fermionic Green's function including point 
particle path integrals, energy momentum integrals, and string modifications to the 
Green's function. We also investigate the relation of the Wick rotated Green's func- 
tion relation to discrete field theory. In section 9 we discuss the main conclusions 
of the paper. In Appendix A we discuss the permutation of the position, time, mo- 
mentum and energy and how this effects the physical picture of the zeta function. 
In appendix B we discuss the one particle partition function in more detail. In 
appendix C we discuss a similar treatment for the scalar particle. 



2 Review of Dirichlet functions 



In this section we briefly review various Dirichlet functions used in this paper. For 
more details the reader can consult texts such as [5], [6], [7], [8] and [9]. Dirichlet 

oo ] 

functions are defined through series expansions of the form f(t,a)= E a n tt • 

n=i n u+u 

Some prominent examples of Dirichlet functions are the Riemann zeta function, 
Dirichlet eta function, lambda function and beta function defined by: 



({it + a) 


oo 

= E n~( lt+a > : 

n=l 


a > 1 




rj(it + a) 


oo 

= E n~t lt+a >(- 
n=l 


_l)n+l 


= _(2 1 -( i *+ CT ) - l)((it + 


X(it + a) 


= g (2n+l)- 

n=0 


-(it+a) _ 


= (i _ 2-( it+ ^)((it + a) 


P(it + a) 


oo 

= E (2n+l)" 

n=0 




-l) n 



In terms of these functions the Riemann hypothesis is equivalent to r/(it + a) ^ 

for — < a < 1 . Dirichlet functions obey a functional relation that relates its value 
2 

at s = it + a to its value at 1 — s = —it + (1 — a) . The extension of the Riemann 
hypothesis from - < a < 1 to the whole critical strip < a < 1 is straight forward 
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using the functional relation for eta: 

r ? ( S )=2vr- 1 sin(^)r(l- S )r ? (l-,) (2.2) 
or in terms of the functional equation for the Riemann zeta function: 

7T S 

C( s ) = 2%- 1 sin(-)r(i- s )C(i- s ) 

C(l-s) = 2(27r)-cos( T )r( a )C( a ) (2-3) 
n-^T(s/2)C(s) = n-fi-ytrUl - s)/2)C(l - s) 

In addition to the series representation for Dirichlet functions one has integral 
representations. For example the eta and Riemann zeta function are given by: 



n( 8 ) = / d£ — ^ -e s£ = / -w 3 - 1 : < a < 1 

,w r(s) 7 exp(e- s ) + l r(s) 7 expw + 1 

o 

(2.4) 



and 



oo oo 

C(s) = -i- / d£ — ^ e s£ = -i- / u^ 1 : a > 1 (2.5) 

w r(s) 7 exp(e- £ )-l T(s)7 expw-1 y J 

— oo 

The Dirichlet functions can also be represented as an integral over theta functions. 
For example the Dirichlet eta function can be written as: 

S oo 

^H^tt^-L- I ±L T s/* { e 4 (o\iT) + e 2 (o\iT)-e 3 (o\iT)) (2.6) 

and the Riemann zeta function can be expressed as: 

S 00 

C(s) = 7r2J^f ^r s / 2 (e 3 (0\iT) -l):a>2 (2.7) 

Several other functions can be defined which place these Dirichlet functions in a 
more generalized context. The Riemann zeta function can be seen as a special case 
of the Hurwitz zeta function defined by: 

oo 

Q H (it + o,a) = Y,{n + a)- { i t+a) (2.8) 

n=0 

which is in turn a special case of the Lerch zeta function defined by: 

oo 

<f>(x, it + a,a) = J2(n + «)-(**+-) e 2 ™* (2.9) 

n=0 
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The Lerch zeta function also has the integral representation for a > : 

oo 

<j>(x,it + a,a) = —r r / dEe~ itE ~ aE — = — ; =r- 

y ' T(it + a) J exp((a - l)e~ E + 2mx) - exp(ae~ E ) 

— oo 

oo 

= / dwW S ~ l r 

l(s) J 1 — exp(— w + 2ttix) 

o 

(2.10) 

and in terms of theta functions is expressed as: 

T s,2 * 



(x, s, a) + e~ 2nix <f>{-x, s, 1 - a) = f drr 3 / 2 ^ 2 ^ e~* (™+«) 2 V™* 

r ( S / 2 ) J neZ 

(2.11) 



The functional relation of the Lerch zeta function is: 

^2-ninx 

4>(x, l - s,a) =J2 



n in + a) 1 s 
1 1 

2m(-S - ax) e -2mna 2m(-S - ax) e 2nina 

(27r)- s r( S )e l 4 E 7 — ; — + (2vr)- s r(s)e l 4 £ . 



n (n + x) s _ n (n + l-x) s 

27ri(-S — OX) 2ni(-S — «l) 

= (2^)- s r(s)e 4 >(-q,s,x) + (27r)- s r(s)e 4 ' <j>(a,8,l - x) 

(2.12) 

Another useful function is given by Laurincikas as: 

oo 

v ( x ,s) = Y / e 27Tixn n- s (2.13) 
n=l 

for a > 1 with integral representations: 



oo 



s/2 \ " — 7TT(n) 2 27rinir\ 





■oo 



n>0 



w(s,i) = 4r / d£e~ aB ; — ; — gr- (2.14) 

^ ' ; T(s) J exp(2vrix - exp (e~ E ) 1 j 



""(s) j exp(27rix) — exp(e E ) 

— oo 

oo 

1 ' dwK" 1 7 — \ ) 



T(s) J exp(2nix) — exp(w) ' 

This function is given by the Lerch function for a = 1 as ip(x,s) = (f>(x, s, 1). 
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In terms of the Lerch Zeta function the Dirichlet functions are given by: 



4>{x = 0, it + a, a = 1) = ((it + a) 
4>(x = - , it + a, a = 1) = — i](it + a) 

2 1 r*M ^ ( 2 - 15 ) 

<f>(x = 0,it + a,a = -) = 2- {lt+a ^ \{it + a) 

<t>{x = i it + a, a = -) = 2- (it+a ^{it + a) 
and in terms of the Lerch function the Riemann hypothesis is 

0(it + a, x = -, a = 1) = -ri(it + a) ^ (2.16) 

for - < cr < 1. £(s) and /3(s) are examples of Dirichlet L functions and are denoted 
by L + \ and L + 4 respectively. Other Dirichlet L functions can be defined from linear 
combinations of the Lerch function at fractional values of a. Clearly the Lerch zeta 
is a key mathematical object. In the next three sections we describe how to represent 
this zeta function and its parameters as a transition amplitude of a 2D fermionic 
quantum field theory. 



3 Propagator, transition amplitudes and Green's func- 
tion associated with the Lerch zeta function 

In the one particle description the propagator, transition amplitude or Green's func- 
tion of a quantum system is represented by: 

A(l 1) = G(x,t;x',t') = (x,t\ x',t') = (x' \U(t-t')\ x) = (x' | e -^(*-*0| x ) 

(3.1) 

Here E is the one particle energy operator. In a second quantized description one 
represents this Green's function as the two point function: 

G(x, t; x', t') = (0\T^(x, t) i)(x\ t')\ 0) (3.2) 

where 4>(x,t) is a second quantized field. The wave function is promoted to an 
operator and | 0) is the vacuum state. In this paper we will always choose tf < t 
so that the expression will be implicitly time ordered. The one particle and second 
quantized descriptions are related because the one particle state is given by | x, t) = 
ip'(x,t) \ |0). For a free particle, the transition amplitude is completely determined 
by the dispersion relation between energy and momentum as well as the spatial 
boundary conditions and whether it is a boson or fermion. 

We consider four cases of dispersion relations and derive the associated fermionic 
Green's functions. These are: 
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(1) E = p , a right moving particle moving at the speed of light. 

(2) E = p , a right moving particle with x taking values on a spatial circle S* 1 of 
circumference L = lixR . 

(3) E = — log(ap + 1) , a right moving particle with logarithmic dispersion and 

(4) E = - log(ap + 1) , a right moving particle with spatial coordinate x taking 
values on a spatial circle S 1 of circumference L. 

The first two cases are standard and are included for pedagogical purposes 
Case (1): E = p 

The dispersion relation E = p with p > represents the usual case of a right moving 
massless particle on space-time R x R. Classically we denote E as the one particle 
energy, p the momentum, v the velocity, x position and t time. The basic classical 
relations are given by the formulas: 

Energy : E = cp 

E 

Momentum : p = — 

Velocity : v = ^— = c (3.3) 
op 

Position : x = ct 

Time : t = —x 

c 

The dispersion relation is given by the first two lines in (3.3) and depicted in figure 
1. 



Figure 1: Dispersion relation for Case (1) E = p or E = — p 

For the most part we shall choose units where c = 1 and /i = 1 however some- 
times we will include them for illustration. Quantum mechanically the dispersion 
relation becomes the operator equation E = p . The momentum operator then 
has momentum eigenstates denoted by p \p) = p \p). In the position basis we have 
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spatial wave functions: ip p (x) = (x\ p) = e ipx . The time dependent wave function 
is given by: ip p (x,t) = (x\e~ lHt \p) = e ~ lE P t e t P x = e -tpt e ipx an j transition 
amplitude is expressed as: 

Ai = (x',t'\ x,t) = (x 1 \U(t'-t)\ x) = (x' \e- iH ^'-^\ x) = Y,e~ iEp(t ~ t ' ) ip{x)^*(x') 

v 

(3.4) 

Using the specific form of the wave function in the position basis we have: 



Ai = ^ 9l (x - x',s) = G 1 (x,t;x',t') =J2e- lEk{t '- t) Hx)r(x' 



k 

oo 

— / dpe -ip(t-t')+i P (x-x') = 1 i 

2vr J 2tt(x- x') -it- t') 





(3.5) 



So that: 



g 1 (x,s) = 2irG 1 (x,t;0,0)= . (3.6) 

x + is/ a 



where s = it /a. 

To see that this is the correct form of the Green's function note that Wick 
rotating to Euclidean space G(x, i, 0, 0) = d x _nK(x, i, 0, 0) = 2d x _ i i— log(a; 2 + 

t 2 ) = — a standard result in 2d field theory. Here K(x,t;0, 0) is the scalar 

two point Green's function and G(x, t; 0, 0) is the two point chiral fermionic Green's 
function. Wick rotating back to real space we obtain (3.6). Note the usual notation 
for G and K is S and D , we don't use this notation as it sometimes conflicts with 
notation for the action and Dirac operator. 

To obtain the Green's function in the second quantized field theory formalism 
one can write the Dirac equation for a chiral fermion in 1+1 dimensions as: 

(i±-)(*).o (3?) 

We use the representation of Gamma matrices given by: 

'-(: oM: o—Mi-i) (38) 

The Dirac equation reduces to (dt + d x )ip = with solutions i^(x,t) = e -*P*e*P x . 
These wave functions are normalized by: 

(ip,x) = J dxtfloX (3-9) 
In the second quantization form of the propagator the Dirac field becomes an oper- 
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ator with mode expansion: 



$(x,t)= J dpb p e~ itp+ixp (3.10) 



where the operators b p obey the anticommutation relations {b p ,b p >} = 5{p + p') 
The second quantized form of the Green's function is: 



(3.11) 



Gi(x,t,0,0) = (0|7ty(a;,t)V>(0,0)| 0) = £ (0\b p e- itp+ixp b_ p ,\ 0} 

p,p'=0 

OO OO flrf. \ j 

_ ^2 g— itp+ixp _ J ^£_ e ~ttp+txp _ j; 

P =o o 27r 27r - *) 

in agreement with (3.5). 

The Weyl Dirac equation follows from the action: 

1 = J dtdx^dt + (3.12) 
and the Hamiltonian derived from (3.12) is given by H = J dxipj 1 d x i/j . In Fourier 

oo 

space the Hamiltonian is written H = J dppb_ p b p . . 

o 

Case (2) E — —— 

Lj 

This is identical to case (1) except the spatial direction takes its values on a circle 
and space-time is a 1+1 dimensional cylinder R x S l . As is well known on such 
a space-time the momentum is quantized as the wave function obeys the periodic 
boundary condition ip(x + L) = ip(x) . Then the basic equations are: 

c2irn 

Energy : E = cp = — — — 

■ dE 

Velocity : v = — — = c 

dp (3 13) 

Position : x = ct y ' 

Time : t = - 

c 



Classically all particles are right moving and travel with the velocity of light. If 
a set of particles leave position and time they all arrive at position x at the 
same time regardless of energy. In the quantum case the above dispersion relation 
becomes the operator equation E = p. The momentum eigenstates are denoted by 
p\p) = p In the position basis we have: (x\ p) = ^ e l?I . Now xeS 1 so that 
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momentum is quantized p = —— with time independent staes 



2irn 

iP n (x) = (x\p = ^\ =-^e~ X (3.14) 



L I LV2 
which satisfies the periodic boundary condition 

ipn(x + L) = ip(x) (3.15) 
The transition amplitude or Green's function is then given by: 

A 2 = ^92{x - x', s) = G 2 (x, t; x' , t>) = £ e^*'-^^)^') 
2irn , ,. 27m . 

™=° e L >( (a: _ x ') _(*_*/)) (3-16) 

I — e L 

1 1 



" 27ri? ^-^-(t-f)) 
1 — e -ft 

so that the Green's function is written: 

<7 2 (x,s) = 27rG 2 (x,t;0,0) = i (3.17) 

i-(x + is/a) 
1 — e -K 

with s = it I a . Note as L — > oo the Green's function reduces to that of Case (1) 

Gi(M;0,0) = i ^7- 

2ir x — t 

In the second quantized formalism the chiral Dirac equation is given by: 

1 ; 2 I ) (3.18) 



with solutions given by: 

2irn . 2-kti 

1 -it— h IX—— 

4>(x,t) = -^e L L (3.19) 

These solutions obey the periodic boundary condition ip(x) = ip{x + L) . The mode 
expansion is: 

oo _ it ^L + ix ^ 
4>(x,t) = Y / b n e L L (3.20) 

n=0 
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where the operators obey {&„,&_„/} = 5 nn i . The Hamiltonian is 

H = — J2 n b-nb n (3.21) 



The Green's function in the second quantized form is given by: 



27m lirn 



1 ~ , , -it—— + ix- 



G 2 (x,t, 0,0) = (0\Tip(x,t)ip(0,0)\ 0) = y V (0|6 n e L L &_ m | 0) 

m,n=U 

^ 27m 27rn 
1 ~ -a—;- + IX-p- 1 1 

™=0 -«*-r- + — 

1 -e L L 



(3.22) 



in agreement with (3.16). 



Case (3): E = - \og(ap + 1) 
a 



In case (3) the energy varies logarithmically with the momentum. This is a non- 
standard dispersion relation. However it reduces to the standard case (1) E = p in 
the limit a — > . The relation of (E vs p) is plotted in figure 2. 




Figure 2: E vs p for Case (3) logarithmic dispersion. 

We shall mainly be concerned with the upper right quadrant associated with a 
right moving particle with positive energy. We can think of a as a new length scale 
that modifies physics at short distances. 
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The classical theory is described by the equations: 



Energy : E = — \og(ap + 1) 
a 1 

Momentum : p = —{e aE — 1) 

dE a c 

Velocity : v = -—- = = ce~ aE 

op ap + 1 



(3.23) 



Position : 
Time : t 



x = 



-t = e~ aE ct 



ap+ 1 

a P +1 x = e aE* 
C C 



These equations yield the nonintuitive relation that the more energy or momentum 
one puts into a particle the slower it goes. That velocity is inversely related to 
momentum and vanishes exponentially with energy. So if a set of particles are 
released at time and space it will take exponentially more time to arrive at 
point x depending on the energy. Nevertheless as long as one stays at energy and 
momentum far less than - the usual picture of case (1) emerges. In figure 3 and 4 
we plot the velocity as a function of momentum and energy. 




Figure 3: Velocity as a function of momentum for Case (3) v = 



ap + 1 



0.8 



. 6 



. 4 



0.2 




Figure 4: Velocity as a function of energy for Case (3) v = ce 



One way to motivate the nonstandard dispersion relation is the following. Con- 
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sider first the Lagrangian L that is sometimes used as a starting point for the 
description of a fermion on a Euclidean semilattice aZ x R. 

L = ip(x,t)(-(ip(x,i+a) -i/j(x,i- a)) + id x ip(x,t)) (3.24) 

We can write this using the exponential of the derivative operator as: 

L = ij(x,t){^(e a9 ^(x,t) - e~ adi eij(x,t)) + id x ^{x,t)) (3.25) 

Now Wick rotating back to real time using t = —it we have formally: 

L = ^(x,t){^(e iadt iP(x,t) - e- ia9t ei;(x,t))+id x ^(x,t)) (3.26) 

After Fourier transforming the fields we obtain the dispersion relation: 

i sinh(a£) - p = (3.27) 

or inverting the relation E = — sinh~ 1 (ap) = — log(ap + \Jl + (ap) 2 ). If instead of 
(3.24) if one started with 

L = i/j(x,t)(-(i/j(x,t + a) - tp(x,t)) +id x Mx,t)) (3.28) 
a 

we can again write this using the exponential of the derivative operator as: 

L = il>{x,t){\e ad ^{x,t) - 4>(x,t)) +id x tp(x,t)) (3.29) 

Now again Wick rotating back to real time using i = —it we have formally: 

L = i;(x,t){^(e ia9 ^(x,t) - ip{x,t)) +id x 4>(x,t)) (3.30) 

Fourier transforming this Lagrangian leads to the dispersion relation 

^(e aE -l)-p = (3.31) 

or inverting the relationship E = — log(ap + 1) as in (3.23). We shall return to a 
euclidean quantum field description in section 8. For now we will continue with a 
first quantized description based on (3.1). 

In the first quantized description the spatial wave functions are the same as in 
Case (1) (x\ p) = e tpx . The time dependence is modified however and from id t ijj = 
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E p ip we have: 

.„ . . -it- logfap + 1) . i— ■ 

$ p (x,t) = e - lE ^e ipx = e a &y F ' e ipx = (ap+ 1) a e ipx (3.32) 

Then using the general formula for the transition amplitude or propagator: 

A 3 = G 3 (x, t; x', t') = £ e-^'-V^x)^ * (x') 
k 

1 °f -i-log(ap+l)(t-t') + ip(x-x') 



= ^Jd P e a 
o 

oo i 

= 2vr / dp ^ aP + ^ 
o 

oo 1 . 

= —-e-^-*')/ 11 f duu~a ~ ' e i{u)i,x-x')/a 
2vra J 



1 *^ e ip(x-x') 



1 

11 ■/ f\ I , i-(t-t')-l °f -i-(t - 1') + 1 - 1 

^7T (X ,/ 



e 

-i((x—x')/a 



= ^i e -i(x-x')/a u x _ a , / y a ^ a ( t *) 1 r (_j(( x _ x ')/ a> _jl ( t _ t ') + 1) 
27r a a 

(3.33) 

where 

oo 

r(M)= J dyy^e^ (3.34) 

and we have used the change of variables u = ap+ 1 and y = —i(x — x')u. Defining 

it 

g 3 (x, s) = 2ttG 3 (x, t; 0, 0) for s = — we express write the Green's function as: 

g 3 (x, s) = 2ttG 3 (x, t; 0, 0) = I e -«/«(I^)-i+«r(^, i _ s ) (3.35) 

Case (3'): E — -log(ap + l) using (E,p) representation 

a 

This is a variant on the above. In this case the Green's function is given by the 
(E,p) representation: 

^-g 3 <(x,s) = G 3 ,(x,t- ia; 0,0) = / dEdp- e -iE{t-ia) jpx 

(3.36) 
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Performing the integral over p by contour integration we have: 



g 3 , (x, s) = 2vrG 3 , (x, t + ia; 0,0)=/ dEe~ iE{ - t+ia ^e i{ - eaE -^ x / a 
= I f dWW~ s e i{w - 1)x/a 



(3.37) 



a J 



where we have defined W = e . Now setting u = —iWx/a this becomes: 




duu s e 



u 



(3.38) 



using s = it we write this as: 




-lg-ia/op^ _ 



(3.39) 



We discuss the (E,p) representation and other representations of the Green's func- 
tion in more detail in section 8. 

Second quantized Case (3) 

In the second quantized description of the Green's function one solves the Dirac 
equation associated to the dispersion relation (3.31) : 




(3.40) 



which reduces to the equation: 



(l( e ^o _ 1)+dx) ^ = 



(3.41) 



Alternatively one can consider an equation of the form: 



(id t + \og(aid x + l)tp = 



(3.42) 



In either case we have solutions of the form: 




(3.43) 
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To second quantize one forms a superposition of these solutions and promotes the 
coefficients of the expansion to operators as: 

00 1 . . 00 

1 r ~ -it— m(ap + 1) + ixp 1 r ... 
^(x,t) = -= j dpb p e a = 7m J d PM a P+ 1 ) /a e txp 



—00 —00 



(3.44) 

As in (3.10) b p> Q are annihilation operators, 6 p <o are creation operators so that 
operating on the vacuum state b p \ 0) = 0(—p)\ 0) and (0| b p = (0| 0(p). The Green's 
function is then: 

2vrG 3 (x,t,0,0) =2vr (0\Tip(x, t)ip(0, 0)| 0) 

00 00 
= J dpdp' (0| 6 p b p /| °) ( a P + l) _i */ a e^ = / dp(ap + i)-«/a e «P (3.45) 
-00 ' ' 



in agreement with (3.33). 

The non standard Dirac equation (3.40) follows from the action: 

/= f dtdx$(~f°d t + 7 1 H-log(l + aid x ))ip (3.46) 
J a 

The differential operator 

D x = -i-log(l + aid x ) (3.47) 

is defined by the series expansion: 

-. 00 

- log(l + aid x ) = ^"V-^-lf- 1 ^ (3.48) 

n = l 



-I 

a 



An alternative action that yields the same solutions to the equation of motion is: 
I=f dtdxi;(~f — (e aidt - 1) + ^d x )^ (3.49) 

J CI 

where the differential operator 

D t = ^(e ai9t ~ 1) (3.50) 
is defined through the series expansion: 

00 -, 

— (e aidt ~l) = J2 a"" 1 *" -1 — ^ ( 3 - 51 ) 

a n=l ' 
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Case (4): E = - \og(a^ + 1) 

(X Li 

This is the case that is relevant to the Riemann hypothesis. As in case (2) periodicity 
for the x field means that it takes its values on S 1 of radius L and this implies 
p n = —— . The dispersion relation is depicted in figure 5. 



5 10 15 20 25 



Figure 5: E vs p for Case (4) Logarithmic dispersion with discrete momentum. We 
have assumed R = a. 

Again we are mainly concerned with the upper right quadrant. The spatial 

2-Kn 
1 i——x 

wave functions are i/) n {x) = (x\ p n ) = — j /9 e L and time dependence of the 



wave functions is given by: 



1 



L 1 / 2 



-iE„t. 



2irn 



[n 



Mx,t) = (x,t\ Pn ) - (l)1/2 . _ (L)1/2V . L 

simplicity we set x' = t' = . Then the propagator is written as: 



2im 

2tT ti i - x 



At = G 4 (x, t; 0, 0) = (x\e- itH \ x = 0) = £ e 

X n=0 
= g (a— + 1)-Wa e i2nnx/L 

27F L n=0 L 

1 2vr / L \ l */ a — / r, \ - %t l a 



-iE n t 



V'n^)C(O) 



27r L \27ra 



71=0 



L 

2ira 



i2irnx/L 



1 2tt / L Y t/a U L ) 
2ir L V 27ra / L ' a ' 27ra 

J_If^)«/aA f _E_ !! ^ 
2vr i?V ^ l 2vri?'a'a J 



(3.52) 



it 

In terms of s = — the Green's function becomes: 
a 



94 (x,s) = 27rG(x,t;0,0) = h-)^ 

it a 



, x R. 
^2^R ,S, ^> 



(3.53) 



So that we have obtained a physical representation of the Lerch zeta function 
(f)(x, it + a, a) where we can interpret the parameter t as time, x as space and a 
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proportional to radius of the universe ——= R . We discuss the interpretation of 
the parameter a in section 4. 

In the second quantized description of the Green's function one solves the Dirac 
equation associated to the dispersion relation (3.31) : 

( 7 oI( e ^o _ 1} + ^ (i±73) ( J ) V = (3.54) 

which reduces to the equation: 

(^(e atdo -1)+ 0*^ = (3.55) 

with periodic boundary conditions ^(x + L,t) = tp(x,t) and with solutions of the 
form: 

27rn 27rn 

j -iln(a h lit i X 

V>M) = ^ e L (3.56) 

To second quantize one forms a superposition of these solutions and promotes the 
coefficients of the expansion to operators as: 

^ 1 , 2im . 27rn 

? 1 ^ £ log a— + 1 ) + zx— 

j n= ~°° 1 2™ (3-57) 

2^ MO-T + 1 ) ae L 



(L)V2 



n=— oo 



6„>o are annihilation operators, 6 n< o are creation operators so that operating on 
the vacuum state b n \ 0) = 6(—n)\ 0) and (0| b n = (0| 9(n). The Green's function is 
then: 

G 4 (M,0,0) = (0\Tip(x,t)i>(0,0)\ 0) 



27T77. 

i-| £ ( 0|M,,|0)(aHp + l)-<./« e "— 

n,n'=— oo 
oo 27rn 

n=0 



(3.58) 



In agreement with (3.52). 

Interpretation of t for Dirichlet series 

Using the formula (3.52) and (3.53) various Dirchlet series can be given a physical 
interpretation in terms of transition amplitudes. For example the Dirichlet eta 
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function is given by: 

G(x = |, t; 0, 0) R=a = /x = | , t\x = 0, t = 0) R=a 

i oo 1 (3.59) 

n=0 

when a = = 1 and Riemann zeta function given by ((it/a) = (2 1 ~*'/ a — 

2ira 

l)r]E>(it/a). Here t is given the interpretation of time and the one dimensional 
volume obeys the condition L = 2ira or R = a. 

The rest of the Dirichlet functions discussed in section 2 are given in terms of 
the Green's function through: 

G 4 (x = 0, t; 0, 0) R=a = = 0, -, a = 1) = ^C(it/a) 

1 x 1 it 1 
G A (x = ira, t; 0, 0)^ = - — 4>(- = -, -, a = 1) = -- — r}(it/a) 

lira L 2 a 2ira 

G 4 (x = 0, t; 0,0) a = — <Ky = 0, tf, a = -) = _2-«/°A(i*/a) ( 3 - 60 ) 
_r= — Z7ra 1/ 2 27ra 

G 4 (x = vra, i; 0, 0) a = = Lit, a = \) = ^-2- a [3(it / a) 

r= — 2ira L 2 2 2na 

2 

Introducing a through t — > t + io we relate the Dirichlet functions to the Green's 
function at analytically continued time. 

Generalization with constant vector potential 

In the presence of a gauge potential A\ one can use the substitution p — ► p + qlAi to 
see the effect on the Green's function. For a spatial S 1 a constant gauge potential 
cannot be gauged away because of the Aharonov Bohm effect on a nonsimply con- 
nected space. The dispersion relation (3.31) is then modified to E n = - \og(ap n + 1) 

a 

27rn 2ir 
with p n = — h qA\ so that p G — Z + qA\. The transition amplitude or Green s 

lj Li 

function is modified to : 

1 9-7T °° 

(x,t \0,0) = ——J2 e~ iEnt+iPnX 



2vr L „ 



1 2tt f L _y*/ a gi;C(?Ai ~ / i LqAi \ lt/a j2Knx/L (3 61 ) 



2vr L V2vra 



n=0 



= i?!fiy 1/fl e^^ !* - 1 LqAl ) 
2ir L \2ita) V L' a ' 2vra 2vr 7 

In this case the boundary condition on the wave function becomes ip(x + L) = 
e iqAlL ijj(x). Note that in two dimensions the charge q has dimensions of [q] = 
(length) -1 while the vector potential A\ is dimensionless. In the presence of constant 
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A\ the value of the parameter a in the Lerch zeta function is a = 1 — - = 

2ira Z7T 

L(— + — ). 

An especially useful form of the dispersion relation in the presence of a constant 
gauge potential is written as: 

R 

Where we have defined A = RqA\ and A' = h A. Two important cases are (i) 

a 

R 1 1 R 

A' = 0, A = and (ii) A' = —,A= which we shall use in section 4. The 

a K 1 2' 2 a 

dispersion relation can be used in deriving the (E,p) form of the Green's function 

given by: 



(n + A') e aE 



Gr(x ,,0,0)U,« = ^ Ejf-^-j^y*- (3.03) 



n=— oo 



i? 7 a 2 
Introducing the proper time parameter C we have: 

G 4 '(x,t; 0,0)1^' 

™ , „ (n + A'V 9 e 2ajB , n 

_J_ v / d c^(<!i±i2 + f!! )e - c «Sr J) » e 'iiV« 

2tt.R ^ 7 2ir y R a ' 

n=—oo 

(3.64) 

Now defining = e aE we write this as: 
G A ,(x,t;0,0)\R, A , 

_ J- V /,i C j5L ( (!l±^ + ^) e - C,< HH' "^'e'iiV-'- 
2vri? ^ 7 2vraTy v i? a 

n=— oo 

(3.65) 

Further defining u = —Ca~ 2 W 2 the Green's function becomes: 

G 4 '(M;0,o)|fl,^ 

= — E / d " (n + ^ ) e" C((i! ^ i])2 e-»e t i X (-n)-^(Ca- 2 )^ 

rt=— oo 



+-L £ / ^-^^^e" C((_ ^ } e-«e^ X (-n)-^(Ca- 2 )^ 
27T.R f-^ y 4vrau G 1 / 2 v 1 y ' 



n=— oo 



(3.66) 

Performing the integral over it and C and expressing the sum over n in terms of the 
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Lerch zeta function we have: 



27riG 4 /(x,t; 0,0)1^' 

it x 

/? — -)- 1 1 T — 2"7n T 

= ( a ] a W * + A>) ~ 6 27lR 0( ~ 2^' S + hl ~ A ' )] C ° t(7rs) 



(3.67) 



4 Interpretation of the parameter a 

We have found an interpretation of the parameters t, x, a of the Lerch Zeta function 
in terms of a fermionic Green's function. What about the parameter a ? Certainly 
this is a crucial parameter for the Dirichlet functions as the Riemann hypothesis is 
given by (f)(x = -, it + a, a = 1) 7^ for - < a < 1 and the functional equation 
relates a to 1 — a . 

One way to interpret a is in terms of scattering amplitudes. After all the two 
point fermionic Green's function is a special example of a (1 — > 1) scattering am- 
plitude. Scattering amplitudes because of their analytic nature can be continued in 
the complex plane and even to other signatures of space-time as have been remarked 
in [10]. 

Another way to interpret a is in terms of a one particle partition function. This 
is related to the Green's function by: 

z ne(/3 = t) = J dxG(x,T;x,0) 

n=o a ^ a 

where r is a Wick rotated time. 

One can also study the partition function in a Harmonic oscillator formalism 
nonstandard logarithmic oscillator with energy E = - log(aw(n + -) + 1) . In this 

CI £ 

case the partition function is: 

^osc(/3 = t) = J dxG osc (x, r; x, 0) 

= £ e-P E " = {au)-P/ a g (n + — + \)~Pl a (42) 

n=0 n=0 aL0 2 

Again because of the analytic properties of the partition function one can study the 
function in the complex plane and relate its zeros to phase transitions of the theory 
in the sense of Yang and Lee. 

However one does lose some intuition when introduces a through analytic con- 
tinuation. After all when some asks what time has elapsed one does not respond 
10hours + i5minutes, or if someone asks what the temperature is one usually doesn't 
come back with (.5 + H4.3)°K. We shall return to the statistical mechanics associ- 
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ated with the Dirichlet series in section 5, and we shall discuss the relation of the 
Harmonic oscillator formalism to the fermionic description of the 2D string matrix 
model in section 6. In this section we will investigate some other physical interpre- 
tations of the parameter a to obtain some insight into the effect of the parameter 
without using direct analytic continuation. 

4.1 a as a mixing parameter 

To begin with consider the first quantized description of the transition amplitude. 
The superposition of momentum space is used to define the states: 

oo 27Tn 

b) = E e ' L X \Pn) ( 4 - 3 ) 
n=0 

Consider a similar set of states parameterized by a parameter a through: 

oo 

k) = E e ~ aEn I Pn) = \ X = °) ( 4 - 4 ) 

n=0 

oo 

and | a) has been normalized so that {a\ a) = z{2a) where z(2a) = J2 e~ 2u ■ The 

n=0 

probability of measuring a momentum value p n from an initially prepared state | a) 
I (p n \ a) I 2 e ~ 2aE « 

is given by P n = — j^-. — ; — = — -. — — . These probabilities satisfy < P n < 1 and 
(<7 | a) z(2a) 

oo oo 

TPn = —7- — r y e~ 2 ° E - = 1 (4.5) 

z(2a) ^ V ; 

n=0 v ' n=0 

The first quantized description of the transition amplitude from: 

00 

A H = (x,t\ a) = (x\ e~ im e~ aH \x' = 0) = ^ e~ itE "~ aE " e ixpn (4.6) 

n=0 

In the second quantized description we use the correspondence: 



\x, t)^fP(x,t)\ 0) 

e- aH \x, t) = \x,t- ia) -> ^(x, t - ia)\ 0) 



(4.7) 



Then the Green's function becomes: 



G a (x, t; x', t') = G(x, t; x',t' + ia) = (0\ip(x, t)ip(x',t' + ia)\ 0) 

= (o\ip(x,t)e aH 'ip(x',t')e- aH \ 0) = (o\e- aH i>(x, t)e aH ip(x', t')\ 0) (4.8) 

= <0|^(x, t - ia)ip(x',t')\ 0) = G(x, t - ia; x', t') 
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oo 

where H = J2 E n b- n b n is the second quantized Hamiltonian. So in the second 

n=l 

quantized form the introduction of the mixing parameter a is simply realized by 
evaluating the Green's function at complex time t — ia . Note the mixing param- 
eter dependent Green's function is not the same as the finite temperature Green's 

oo 

function which is given by: Gp(x, t; x' , t') = J2 (— l) n G(x, t — ni{3; x' , t') although 

n=— oo 

they both involve the complex time. Further discussion of statistical ensembles and 
finite temperature Green's function will be given in section 5. 

Now we specialize to the Cases considered in section 3. First consider case (2) 
and (4) 

Case (2) E = ^ 

For the usual case Hr 2 ) we have E n = p n and p n = — — . We define mixing 

±j 

parameter states by: 



%a / 2) = £ e iPnx e -oE n /2\ pn) 



n=0 

oo 



\x\oj 2) = J2 e^e-^"/ 2 ! p n ) 

n=0 



(4.9) 



These states are normalized through: 

oo -. 

(x , a/2\x, a/2) = (x>, a/2\x>, a/ 2) = z(a) = ^ e"^ = ^ (4.10) 

n=Q 

1 - e L 

Because of the simple time evolution of the momentum states, the time dependence 
of \x,a/ 2) is : 

oo 

\x,a/2, t) = J2 e ip " x - iEnt - aEn/2 \ Pn ) (4.11) 

n=0 

So that the transition amplitude becomes: 

OO -, 

<*' , a/2\x, a/2, t) = £ e*»<?-* ')-^t-aE n = _ _ __ (412) 

n=o i — x — i — t a 

1-e L L L 

The transition probability is then: 

| (x' ,a/2\x,a/2, t) | 2 



P( X ,t^x') a 



(x ,a/2\x,a/ 2) (x' , a/2\x', a/ 2) 

2-ko 



(x>,a/2\x,a/2,t)\ 2 (1 - f L f (413) 



z{a)z(a) 2tt .2vr 2ir 

v ' y ' i—x — i—t — a 

II — e L L L |2 
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This can be further simplified to: 

27r a 

P(x, t - x') a = (l-e L l_ __ (4 m) 

(l-2cos(-^(ar-x' -t))e L a + e L° ) 
The minimum vale of the probability is then given by: 

27T 

Pmin = (1 " 6 2^ ^ = (C0th(^)) 2 > (4.15) 

(i + e "T CJ )2 

This is strictly positive except in the trivial case of a = when there is no mixing. 
Case (4): E = i log(ap„ + 1) 

1 27rn 
For the unusual case we have E n = -log(ap n + 1) and p n = — — . The mixing 

parameter states are defined as: 



2irn 

\x,a/ 2) = g e iPnX e -*E n /2\ pn) = g X {a ^ + ir a/2a| ^ 

n=0 n=0 ^ 

2-7rn . 

\x', a I 2) = E e^-e"^/ 2 ! p n ) = E e L (a— + l)- CT / 2a | p n > 

n=0 n=0 ^ 

These states are normalized by: 



(4.16) 



(x ,a/2\x,a/ 2) = (x' , a/2\x', a/ 2) = z(a) 

= £«-*«•= + = (^r'»&(-,^) (4 - 17> 

n=0 ra=0 ^ ^ 7r O O 

Because of the simple time evolution of the momentum states, the time dependence 
of \x,a/ 2) is: 

oo 

\x,a/2, t) = J2 e iPnX - iEnt ~ aEn/2 \ Pn ) (4.18) 

n=0 

The transition amplitude is: 

oo 

{x ,a/2\x,a/2, t) = E e*-^ - * )~ iE ^-^n 

n=0 

27rn , 

= g e~ (X " X \a^ + l)(- it -^ a = (^)^)/V(— ,— , — ) 
n= o L 27ra L a 27ra 

(4.19) 

where is the Lerch zeta function. Again the transition amplitude is proportional 
to the Lerch zeta function however this time it is evaluated at s = i(t — ia) = it + a. 
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The transition probability is then: 



(x ,a/2\x,a/ 2) (x' , <r/2|x', a/ 2) 



(C»(-,^-)) 2 

a 2-Ka 

Now we can give a physical representation of the Riemann hypothesis by special- 
ly L R 

izing to x — x' = — = ttR and = — = 1. In that case the transition probability 

2 2ira a 



becomes: 



1 it + a _ N , 2 



••1) 

P(x,t^x') = 2- a 



(4.21) 



(Ch(^,1)) 2 

_ |r ? (^)| 2 _ |(l-2M^)/ a)C (!i±^)| 2 

(C(f)) 2 (C(^)) 2 

and the Riemann hypothesis is the statement that: 

|( 1 _ 2 i-(^)/a)£(^±^)|2 
P(x, t ->■ x') = = 2 > o (4.22) 

(C#) 2 

la , L L 

for - < — < 1 and x — x = —, — = a. The physical setup is illustrated in Figure 

2 a 2 2-k 

6. The particle is measured at initial position x = 0. The probability that a particle 

will subsequently be measured at position x = L/2 is nonzero according to the 

Riemann hypothesis. 

x=L/2 




Figure 6: Physical setup for transition probability interpretation of the Riemann 
hypothesis. 
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One can also define the time averaged transition probability 



P < T > = f dwl dn < X>/2|M ' al2) |2 = mm* /*'"^» 2 



T T 

, . At + a , 

itr)( ) 

a 

(4.23) 

The literature on zeta function theory contains a great deal of information on these 
averages for large T as they are easier to study than the Riemann zeros. 

A variant on the above analysis is to consider the transition between states: 

n 

oo 2m — X 
I x) = £ e R I Pn ) 

n=0 (4.24) 

OO OO -n v ' 

I a) = £ e-^| p n ) = E (a- + l)-"/ a | p„> 

n=0 n=0 -K 

The I a) state is normalized by 



OO 

(a\ a) = z(2a) = £ e" 2 ^" 

n=0 

_f (.»+!)-»/. .(VfeA^ 
n= a a a 



(4.25) 



The t evolution of the | x) state is straightforward because of the simple evolution 
of| p n ) by the phase e~ tEnt and we have: 

Tl 

|s,t> = 5> i? (4.26) 

?1=0 

The transition amplitude is then: 

n n 

f) = £ e -° E ™e- itE ™eR X = £ («5 + l) ( ~ CT ^ )/a e^ , , 

n=0 n=o (4.27) 

Taking the magnitude of the amplitude and dividing by the normalization of the 
| a) states we have: 

. , x it + a R.. 2 

\(a\x,t) \ 2 _ \(*\x,t)f 

(a\ a) z(2a) ( H (2a/a, R/a) V ' ; 

One can then use the asymptotic expression for the magnitude of the Lerch zeta 
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function for - < a < 1 and T — > oo [8] 



T 

/ |</)(x,t,a)| 2 cit 



= Ch(2<7, a)T + v ; - Q H {2 - 2a, x)T 2 ' 2a + BT 1 ^ log T + BT ff / 2 
2 — 2cr 

with i? a constant to obtain the asymptotic formula for the time averaged quantity: 



T T I / / it ~\- (J i?, ,n 
1 /■ I (cjIx t> I 2 1 } '-) 

_ / rft 1 - 1 ; 1 = — / dt a a 

T7 (a\a) TJ ( H (2a/a,R/a) 

„ 1 , 1 q(2a-2a)/a,x/27TR) T 1-2^ 

2-— tfrr/a,R/a) l 2vra J (4 . 30) 

-- --1 
S(-) a log(T/a) + S(-)2a 

_j (I (X 



( H (2a/a,R/a) 



Specializing to the case relevant to the Riemann hypothesis we have R = a and 

x = — = ttR = ira so that 
2 



(a\x = ira, t) \ 2 | (a\x = ira, t) 



<cr| a) z(2a) 

\ V C±±^)\ 2 ki - 2 i -^+-)/ a )c(^)i 2 (431) 

C(2a/a) C(2a/a) 
The Riemann hypothesis is equivalent to the statement that the amplitude 

(cj|x, t) (4.32) 

1 cr 

is never zero in the range - < — < 1 . The asymptotic formula of the time averaged 

2 a 
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quantity for R = a ,x = ira ,- < — < 1 , and T — > oo becomes: 

2 a 

T . . . T W ** + ^M2 



T7 a a) Tj 



dt- 



a 



T J (a\a) T J ((2a/a) 
a 

1 (2 2 ' 2 a - l)C((2a - 2a) /a) T 1-2^ 
2 _2a C(2^/a) l 2vra J (4.33) 

-- --1 
B(-) a log(T/o) + J B(-)2a 

_i tt 

C(2a/a) 



4.2 a in Partial sums and neutrino mixing analogy 

In this section we will develop an analogy of the transition probabilities of the 
previous section with the phenomenon of neutrino mixing. Partial sums for the 
Lerch zeta function are truncations of the series representation to N terms and are 
defined by: 

N —it — a x 

X it + <T v — ^ 2irin — 

<My, ,a) = V(n + a) a e L (4.34) 

n=l 

One can formulate a sufficient condition for the Riemann hypothesis in terms of 
these partial sums as: 

1 N 

(j) N (x = -,it + a,a = l) = -r] N (it + a) = J2 {-l) n n~ lt ~ CT / (4.35) 

2 71=1 

for < a < 1 (in this formula we have set a = 1 ). Note that this a stronger 
condition than the usual Riemann hypothesis which only involves N = 00. 

A nice feature of quantum mechanics is that it can be defined for finite dimen- 
sional Hilbert spaces so that one can also represent these partial sums as quantum 
transition amplitudes in a system with a finite number of states. Then quantum 

N 

states are defined as | x,a) = e 2mnx ' L e~ EnCr \ n) with special states given by: 

n=0 

I v) a = |0, a/2) = £ e^"/ 2 | n) 

T n ~N (4.36) 

\v') a = \-,a/2) = Y J {-V n e- aEnl2 \n) 

Z 71=1 

The time dependence of the | v) a is easily determines from the evolution of the 
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momentum basis states 

N 



v, t) a = Y, e - aEn/2 e- iEnt \ n) (4.37) 

n=l 



N 

The normalization is defined by {u\ v) = {v'\ v') = eT aEn = zn{&) So that the 

n=0 

transition probability is: 

N 

, 2 | (— l) n e~ ™*~ n<T | 2 

r, ( l\ I WWi *) I n=l 

P N (y -► l/) = = T — 

N ,n „ ^ ' 6ii > 

z N {2(i) - 2 £ e -^+E„) CO s(t(£ m - £ n ) 

= 

1 27rn 1 L 

Using in the dispersion relation E n = — log (a h 1) = — logfn + 1) for = 

a L a 2ira 

R 

— = 1 we obtain a sufficient formulation of the Riemann hypothesis in terms of 
a 

P(y -► v') = | t) | 2 > (4.39) 

where 



P N [y -► i/ ) - | /..i , A |2| 7-7; » |2 - ^73 > (4.40) 



(f>, t) | 2 = Iryjv^ + a )! 2 
v) | 2 | i/') | 2 4(a) 

To gain further insight into the mixing parameter let us consider truncating the 
momentum sum in (4.34) to two terms so that N = 2. 
First define: 

\ x = ^) = ^(lp = °>-b = i» , N 

2/ f' (4.41) 
\x = 0) = -^(\p = 0) + \p=l)) 

Then introduce the mixing parameter a by defining two states: 

| v ) a = |0, a/2) = £ e —^/ 2 | n ) = | p = 0) + 2~ CT / 2 |p = 1) 

T n= h (4.42) 

I = I v ^/2> = E (-l) n e-^/ 2 | n) = \p = 0)- 2'^ 2 \p = 1) 

The normalizing factor is z(a) = 1 + 2 _CT so we set: 

(4.43) 

where we defined 6 = t&n~ 1 (2~ a / 2 ). The time dependence of the v state is easily 



u) a = cos 9\ p = 0) + sin 9\p = 1) 
= cos 9\ p = 0) - sin 9\p = 1) 



29 



expressed as 

\v,t) = cos0e- iEat \p = 0) + sin9e~ iElt \p =1) 
= cos 6\p= 0) + sm62- u \p= 1) 

So that the transition probability is given by: 

P N =2{v -+v') = \ WW, t) | 2 = |(cos#) 2 - 2~ i '(sin0) 2 | 2 
= (cos 6>) 4 + (sin 6>) 4 - 2(cos 9 sin 9) 2 cos(t log 2) 

= 1 - 2(cosflsin#) 2 (l + cos(tlog2)) (4.45) 
= 1 - (sin20) 2 (cos(tlog v / 2)) 2 



This is always greater than zero except for 6 = — or a = which is outside the 
critical strip < a < 1. The transition probability for N = 2 and cr = .75 is plotted 
as a function of time in Figure 7. It is identical with the two flavor neutrino mixing 
transition probability with AE = log 2. 

Pn=2 (t) N=2 a=.7 5 




t 



2.5 5 7.5 10 12.5 15 



Figure 7: Transition probability for the N = 2 partial sum for mixing parameter 
a = .75. 

The N partial sum transition probability is analogous to the ./V flavor neutrino 
mixing transition probability. In figure 8 and figure 9 we plot the transition prob- 
ability for N = 4 and N = 16. Unlike the N = 2 case we do not have periodicity 
of the probability amplitude. The case ./V = oo corresponds to the Riemann zeta 
function and is considered separately below. 

The physical interpretation of these transition probabilities is as follows. One 

prepares an initial state weighted by energy according to a Boltzmann distribution 

L 

with mixing parameter a . One sets up a detector halfway around the universe at — 
. One then counts particles received and records the time they were received. One 
then plots the number of particles received as a function of time and a . If there is 
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2.5 5 7.5 10 12.5 15 
Figure 8: Transition probability for N = 4 partial sum for mixing parameter a = .75 




Figure 9: Transition probability for N = 16 partial sum for mixing parameter 
<r = .75 
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time when no particles are received then the survival probability drops to zero. If 
L 

the particle isn't at — then the probability it is at another position is non zero. 

We saw this at in the simple two state system above. In that case there were 
two positions x = and x = — (note that because of periodicityx = L is the 
same point as x = ). The probability moves between the two states but the 
transition probability never vanished (except in the trivial case when the initial 
state was pure a = 0) . This is exactly the behavior observed in the oscillation of 
massive neutrinos (or any other two state system that has two different energies and 
evolves from mixed superposition to a single state). Indeed the plot of the transition 
probability is familiar from such experimental neutrino studies. 

The N = oo transition probability discussed in the previous section also has a 

neutrino mixing analogy. In a Kaluza Klein description neutrino mixing involves the 

reduction from 5^4 dimensions. The five dimensional momentum ps = m n = — — 

Li 

acts like a mass in four dimensions through the Dirac equation (ij^d^+i^d^ip = 0. 

n 

oo 2-Ki — X5 

Introducing a Kaluza Klein expansion ip(t,Xi,x^) = E ' t Pm('t^ x i) e ^ the 

m=— oo 

Dirac equation becomes {i^d^+^nin^n = 0. The probability then moves between 
the different mass eigenstates with initial mixing dependent on the interaction with 
the boundary of the five dimensional space. 

In our case the Riemann Zeta function condition a = = 1 means that 

2na 

we are essentially in a Kaluza Klein situation for any reasonably long time scales. 
Essentially the radius of the universe in fixed at a small scale proportional to the 
parameter a . As in Kaluza Klein theory we can effect a reduction from 1 + 1 — ► + 1 
dimensions through the expansion 

oo „ n „ 

^ — , 2m — X\ 

Y>(t,xi)= £ *Pm(t)e L (4.46) 

m=— oo 

with two special states defined by: 

oo 

\v,t)= E e - iEnt e~ aEn l 2 \ n) 

So =1 , ' (4.47) 

W) = E (-l) n e- CT ^/ 2 | n) 

n=l 



So that the transition probability is: 

i <y\v, t) 



2 | E ( — \) n e~ iEnt ~ En(T \ 2 

PI l\ \\" I"' "/ I ?1=0 

[V — > v ) - 



,H-')I 2 IH-)I 2 ^» (448) 

z(2a) - 2 E e~ a( - Em+E ^ cos{t(E m - E n ) 
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and the one particle partition function,^ ) = e n is used to normalize the 

n=0 

states. 

The formula for the transition probability depends only on the dispersion re- 

lation. For case (2) E n = —— we have the + 1 dimensional Dirac equation is 

±j 

(1'KTl\ i —i — - — t 

idt H — —J ipn = with solutions ip n (t) = e~ irTLn = e L> . For case 

i 2im 1 

(4) E n = — log(a h 1) = - log(n + 1) these equations become: 

a L a 



Q(e-*-l) + ?p)^(t) = 



(4.49) 

^ n (t) = e - im «t = + I)"**/- = ( n + l)"«/a 

The mixing in either case comes from a Boltzmann distributed initial state. 

One may wonder where the particle goes so that the probability goes up and 

own as a function of time. As in the Kaluza Klein neutrino case the particle is lost 

to the bulk, in other words there is finite probability it will found at positions other 
L 

than x = or x = — . 

To compare with the Kaluza-Klein neutrino oscillation we note that 



"> = E | a D —a\ n ) O-™) 
n=0 



i R\1 

(n+- 

a 



i? 2 

whereas the Kaluza Klein expansion is given by [11], [12], [13], [14]: 

oo 

I v)kk = E , m = =\ n ) ( 4 - 51 ) 

„=o (n + Rfiv) , 2n 

y( ^ +^) 
i i 

So that the form is the same when m = — , fly = — , l^s = and cj = a. with 

a a 

[iy and /X5 are vector and scalar fermionic mass terms described in [12]. However 
the time variation of this state is very different in our case because of the nontrivial 
logarithmic dispersion relation. 

Using the methods of this section each of the cases 1-4 can be generalized to 
finite a with the transition probabilities where a is a mixing parameter. 
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Case (1) E = p 

The v and v' states are denned by: 

\v,t) = J dpe i7Ta Pe- itp ' ap / 2 \ p) 

\u') = j dpe-°Pl 2 \p) (4 ' 52) 

oo 2 

They are normalized by: (y\ v) = {v'\ v 1 ) = J dpe~ ap = — . The transition proba- 

cr 

bility is given as: 



IH l 2 IH ^) I 2 

I l |2 

|</i(7ra,it + cr)| 2 _ ira-(t-ia) _ cr 2 



(4.53) 



1,2 |i|2 "(vra-t) 2 + a 2 



C7 



where we have used gi(x, s) = computed in section 3. We plot the transition 

x + is 

probability in figure 10 for a = .75a. 




Figure 10: Transition probability for Case (1) E = p. 



ft 

Case (2) E = —with R = a 

R 

The v and v' states are defined by: 



v,t)= £ e - itn / a e- an / 2a \ n) 

ST° , (4.54) 

v >) = g e »irn e -<m/2a| n ) 
n=0 
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oo I 

These states are normalized by: (v\ v) = (z/| v') = Y\ e~ na ' a = r - 

n=0 1 e ' 

the transition probability is: 



. Then 



(i/\u, t) 



(v>\ v>) (v\ v) 



\g2(ira, it + a)\ 2 ' \ _ e i(x-t+ia)/a ' 



1 



1 



1 - e-°l a 

~y gi(x— t+ia)/a ' 



(4.55) 



1 - e~ a / a 
and we have used g2(x,s) 



1 - e- CT / a 
1 



from section 3. We plot the transition 



\ _ P i(x+is)/R 

probability in figure 11 for a = .75a We show in section 8 that case (2) is the 

oo 

periodization of case (1) in the sense that g 2 (x, s) = 9i i x + 27ri?n, s) or 

n=— oo 

oo 

G 2 (x,t,0,ia) = G 1 (x + 2irRn,t;0,ia) (4.56) 

n=— oo 

This explains the multiple peaks in figure 11. 




10 



15 



20 



n 



Figure 11: Transition probability for Case (2) E = —with R = a and a = .75. 

K 



Case (3) E = -log(ap+ 1) 

The v and v' states are defined by: 



v,t)=f dp(ap + l)-i*A»-»/2a| p ) 
!/') = Jdpe inap {ap + l)- CT / 2a | p) 



(4.57) 
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oo I 

The normalization is denned as (u\ v) = v') = J dp(ap + l)~ a > a = 

o a - a 

transition probability is: 



. The 



I WW, t) | 5 



l(^V)l 2 IH I 2 



\g^Ka 2 it + a) 

, 1 ,0 



7T 



o" — a 



(4.58) 



,(j -a (-i- + log(7r))(-it + (l-<7))/a 2 
s je 2 T(—i7r,l — {it + a)\ 



where we have used g3(x,s) = e * x / a (— ix) 1 S T(— ix, 1 — s) for s = (it + o)/a derived 
in section 3. We plot the transition probability in figure 12 for a = .75a. 



P 3 (t) 



E=log(p+l) a=.75 




Figure 12: Transition probability for Case (3) E = - log(ap + 1) for a = .75. 



4.2.1 Case (3') 

This is similar to the above case except g 3 > (x, s) = G 3 /(x, t + ia — ia; 0, 0) is defined 
by the (E, p) integral: 

gy (x, s) = (x,t + ia — ia; 0, 0) 

1 

»l 

-{e aE -l)- p + ie 



, , dEd Pl e -i(t+ia-ia)E e ipx 

^) 2J h e aE_^_ n ,. (4.59) 



i a 

-e- ix ' a (— )" 1+ T(1 - s) 
a a 
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The v and v' states are denned by 

If 1 

\ v i~ t ) = TTZvi / dEd Pi 

( 27r ) J L( e "E_ 1 - 



-i(t+ia-ia/2)E\ 



p) 



{e atj - 1) - p + is 



v>,t>) 



If 1 
——2 / dEdp-, 



-i(t' +ia—icr/2)E ina/at 



P) 



(e atj - 1) -p + ie 



(4.60) 



The transition probability is then given by: 

I WW, t) | 2 



I \ u I u i 
\gy(air,(it + a)/a)\ 2 



(u'\ z/)„| 2 | (u\ v) | 2 

|( e (-i W /2+i Q g( w ))(-it+(i-<T))/oj r ( 1 _ ( it + (7 )/a)\ 2 



(4.61) 



M^)| 2 

In figure 13 we plot the transition probability for mixing parameter a = .75a 

p 3< (t) CT= . 75 
0.51 




Figure 13: Transition probability for Case (3') using (E,p) representation for a 
.75. 
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Case (4) E = - log(o- + 1) 
a H 

For R = a this is the case relevant to the Ricmann hypothesis. The u and u 1 states 
are defined by: 

oo 

I v,t) = £ n- lt l a n- a l 2a \ n) 



n=l 

oo 



v>) = J2 e i7rn n- a / 2a \ n) 

n=0 



(4.62) 
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The states are normalized by (is\ v) = (u'\ v 

') = E n- a ' a = ((a /a) . The transition 

n=l 

probability is given by: 



WW, t) 



\g^a,it + G^ V mtlo)/a)\ 2 _ |(1 - 2 1 -( i *+"))C((it + a)/a)\ 2 (463) 



|C(a/o)P |CWa)| 2 |CWa)| 2 

We plot the transition probability for a = .75 in Figure 14. 

p 4 (t) E=log(n+l) a=.75 

0.4 



. 3 



0.2 



. 1 
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Figure 14: Transion probability for Case (4) E = — log(a h 1) with R = a and 

a R 

a = .75. 

We show in section 8 that case (4) is the periodization of case (3') in the sense 

oo 

that gi(x, s) = J2 9s'( x + 27ri?n, s) or 



G<i(x, t, 0, ia) = ^2 Gy(x + 27rRn,t;0,ia) 



(4.64) 



This explains the multiple peaks in figure 14. Note that unlike the periodization 
from case (1) to case (2) one does not have periodization in t. This because in case 
(3 — 4) the Green's function is not a function of (x — t). The Riemann hypothesis 
in this context is the statement that that the transition probability from the state 
v — ► f ' is nonzero for mixing parameters in the range —a<a<a. 

4.3 a as a SUSY breaking parameter 

In supersymmetric quantum mechanics one considers Hamiltonians of the form [15]: 



(4.65) 
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Here b\,b\ are fermionic creation and annihilation operators. In principle U can 
be an arbitrarily complicated function. Often its form it determined by symmetry 
properties. For example U can be a modular function if y represents a modulus or 
compactified scalar field [16]. If one takes the zeta function as U then the symmetry 
property would be the functional equation together with the Dirichlet expansion 
which uniquely determines the form of the function up to a proportionality constant. 
To be more specific we add to the above the Hamiltonian 

where W = U + iV is an analytic function of x + iy. The total Hamiltonian becomes: 

H = H 1 +H 2 

l ,Pl Pv i ti ,i2n 1, dV.-t, 1, 

= 2<^ + j£ + lWl ) + ^ - 2> + * <*» - 2> 

dll dV 

Now using the Cauchy-Riemann condition — — = — — we have: 

ox oy 

1 r?' V 2 dV 
H = -(«!L+*y- + \W(x + iy)\ 2 ) + —(b\h + b\b 2 - 1) 

2 mi ni2 dy 

Setting W(x + iy) = Q(x + iy) and relabeling x by a we obtain the Hamiltonian: 

H = rEs_ + £y_ + K(<7 + i } | 2) + (6 t 6l + 6 t &2 _ 1} (4 . 69) 

2 mi fn2 dy 
Finally taking nil -toowe have : 

1 T?' flV 

H = -(^ + \((<r + iy)\ 2 ) + -riblh + b\b 2 - 1) (4.70) 
l m 2 ay 

In supersymmetric quantum mechanics the Hamiltonian is related to the super- 
charges by H = Qa . This means supersymmetry will be broken unless the 

a 

ground state has zero energy. Thus supersymmetry will be broken if the function 
W does not intersect zero. In Figure 15 and 16 we plot \W\ 2 = \((iy + -)| 2 and 

|W| 2 = |C(<» + 0)| 2 . 

If the Riemann hypothesis is false then supersymmetry may remain as a sym- 
metry of (4.70) assuming a zero occurs for large y. To understand further how a 
manifests itself as a supersymmetry breaking parameter we need to understand how 
the potential is generated. For example finite temperature effects always break 
supersymmetry because the bosons and fermions react differently to finite tempera- 
ture. Bosons have integer momentum in the imaginary time formalism and fermions 
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(4.67) 



(4.68) 



I w I 2 ct=.5 





Figure 16: \W\ 2 = \((iy + 0)| 2 . The absence of zeros according to the Ricmann 
hypothesis leads to broken supersymmetry. 
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have half integer. This leads to potentials that have no zero for temperature T > . 
In that case the supersymmetry breaking parameter is the temperature T which is 
varied away from zero. In another example finite lattice effects can break supersym- 
metry because the translation symmetry is broken by the lattice and supersymmetry 
is related to space-time translations through the superpoincaire algebra. In that case 
the supersymmetry breaking parameter is the lattice spacing. In our case the zeta 
function has both elements of discreteness from the Wick rotated Dirac equation 
(3.28) and momentum quantization from the periodicity of the spatial direction in 
(3.56). 

4.4 (0 — > 0) String amplitudes and (1 — > 1) amplitudes in Riemann 
dynamics 

(0 — > 0) String amplitudes 

One loop String vacuum amplitudes are defined by two dimensional sigma models 
of the form: 

A(0 0) = J De a z DX M e~ I d 2 «<<^ z xa- z x) (4n) 

The sigma model field X M (z,z) maps the two torus T 2 to a target space of the 
form S 1 x M^ 1 where the circle has radius R . The two dimensional world 
sheet metric is related to the modular parameters of the torus through e^e^dzdz = 

(da^da 1 )^ T * ^ ^° ^ and z = ^ + ia°) , z = ^(a 1 - ia°) . On 

a torus the sigma model reduces to an integral over the modular parameter and its 
complex conjugate so that the path integral becomes: 

A(R) = fdT6(T,f,R)P L (T)P R (f) 

F 

+ I^ya _ ^,Vl - H^)2 ( 4 - 72 ) 
= JdTdf(J2e R «' K R a ) Pl (j)p r (t) 

p n 

where Pl,r(t) = trL,n(e~ TM2 ) is a left or right moving string partition function and 
R is a compactification radius. 

It is knownthat global anomalies can lead to zeros in one loop vacuum amplitudes 
of string models. In [17] models were found that had global anomalies in symmetries 
that were not required for the consistency of the theory so these string models still 
were sensible. For example zeros in one loop amplitudes were found in which the 
integrand of (4.71) had an additional modular symmetry of the form: 



Usual : r — > r + 1; r — > — - 

Additional : r — > — — 
2r 



7 (4.73) 
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The additional symmetry usually arises in the form of a discrete exchange of com- 
pactified spaces. This symmetry together with the modular properties of the inte- 
grand leads to the zeros. 

Also zeros in the derivative of the vacuum energy of string models can also 
be obtained at special values of the compactified radii [18]. These are points of 
enhanced symmetry of the theory where the usual U(l) x U(l) of the compactified 
theory is enhanced to SU{2) or SU(3) and certain states becomes massless. These 
are also the location of thermal divergences if the compactified moduli are inverse 
temperature or chemical potentials in the imaginary time finite temperature string 
theory where the vacuum energy becomes the free energy. 

Other features of (0 — > 0) string amplitudes are the existence of the fundamental 
region: 

F = {t:\t\>1- 1 -<t 1 < 1 -} (4.74) 
and the duality symmetry under the transformation: 

a' 

* = Ji ( 4 -75) 

These are very important properties of string theory which differentiate them from 
traditional point particle theories. 

(1 — > 1) Amplitudes in Riemann dynamics 

We have mainly used the mode expansion representation of the Green's function of 
the one particle to one particle (1 — > 1) amplitude. In section 8 we will derive several 
other representation of the Green's function. One of these is the point particle path 
integral representation: 

x(l)=x,t(l)=t 

A(l -> 1) = G(x,t;0,ia)\ R = J DeoDxoD^D^DPDxDEDte' 1 (4.76) 

x(0)=0,t(0)=i<r 

The fields (eo(r), Xo( T )) define a 0+1 dimensional supergravity associated with (su- 
per)reparametrization point particle action, £ and £ are world line fermion fields 
which eventually become the target space gamma matrices. The point particle ac- 
tion is given by: 

/ = T Jdr(Px -Et- f (P 2 - \(e aE - I) 2 ) 

n la (4.77) 

a a 

The fields (x,t) map the interval r G [n,^] to a cylindrical target space where the 
S 1 factor has the radius R (Recall the string 0^0 amplitude amplitude mapped 
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the two torus onto S 1 x M d 1 target space). The periodicity of the target x direction 

Tl 

ensures that the momentum is quantized as P = — . The Tiechmuller parameter 

K 

of the interval is: C = J eo(r)dT We use the symbol C to denote the Teichmuller 

Tl 

parameter for the point particle to avoid confusion with parameterized time in the 
point particle action. 

The point particle Path integral defining the (1 — > 1) amplitude can be written 
as an integral over C by: 

A{\ -► 1) = (2 s - l)ir- s / 2 T{s/2)aG{x = ira, t; 0, -ia)\ R=a 

= (2 s - l)(2 1 - s - l)vr- s / 2 r( S /2)C(s) = / dcc s / 2 -\e 3 {c) - 9 2 (C) - 4 (C)) 

F 

(4.78) 

The prefactor ir~ s / 2 (2 s — l)T(s/2) is introduced to allow simple transformation asso- 
ciated with the functional equation. We leave the details of the derivation to section 
8. Some features can be recognized however. The theta functions occur because of 
the quantization of momentum associated with the periodicity of the target space 
and the form e~ cp2 in the path integral. The unusual C S//2 factor occurs because 
of the e - E ( tt + a ) factor the change of variable W = e aE in the integral over W. 

The usual modular transformation for a point particle on the interval world line is 
C — > —C associated with the discrete symmetry eo(r) — > — eo(— r) and the exchange 
of the endpoints T\ and ti- As shown in [19] it is this world line modular symmetry 
that leads to the discrete space-time reflection symmetry CPT in the Target space. 
There is an additional modular symmetry C — > — associated with the functional 
relation obeyed by (4.78). One can actually derive the functional equation by tracing 
the transformation properties of the theta functions under C — > — and their effect 
on the (1 — > 1) amplitude. Thus the modular transformations of (4.78) are given 
by: 

Usual : C -> -C 

1 (4.79) 
Additional : C — > — 

The additional modular transformation is directly connected with the Riemann ze- 
ros. This is because the Riemann zeta function is uniquely determined up to a 
constant factor by its functional relation and the fact that it has a Dirichlet expan- 
sion. The functional relation is in turn derived from the action under the additional 
modular transformation. The Riemann zeta function is also uniquely determined as 
a product over the Riemann zeros as: 

«'>= 2 ( .-iTO + ./2) n('-fo- / ' <"») 

where b = log(27r) — 1 — ^7 and p are Riemann zeros. As both the position of the 
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zeros and the functional relation uniquely determine the zeta function they must 
be related. If one interprets s as additional data for the (1 — > 1) amplitude on can 
interpret the Riemann zeros as specific choices of this data that cause the amplitude 
to vanish. In the (0 — ► 0) string amplitude one has choices for the internal radii and 
boundary conditions. In the (1 — > 1) point particle amplitude one has additional 
data namely the initial and final space-time positions and the mixing parameter a. 

The usual modular transformation C — > — C of the point particle action on the 
interval has important implications for the fundamental region F . It allows one to 
choose the fundamental region of the point particle as F = {C : C > 0} = [0, oo][20]. 
For example the energy momentum (E, P) form of the Green's function 

iCtp 2 - —(e aE — l) 2 1 
JdCdPdEe a 2V ' ( P + -( e aE - l))e- lEt+lPx 

F -i a (4.81) 

= / dPdE 1 e -iEt+iPx 

P--(e aE - 1) 

would not be obtained if one chose F = [—00,00]. In that case one would have 
energy momentum delta functions instead of denominators. This would mean that 
point particles would not be able to propagate off shell. We discuss the (E, P) form 
of the Green's function in more detail in section 8. 

a 4 

The additional modular transformation C — > — means that using the theta 
function representation (4.78) we can further restrict the fundamental region as 

F = {C : C > a 2 } = [a 2 , 00} (4.82) 

and we have restored the natural units of the parameter C. The fact that a point 
particle theory can have a fundamental region away from the ultraviolet divergence 
C = is somewhat unexpected as such properties are thought to be exclusive 
properties of string models leading to their finiteness. 



4.5 Duality symmetry and the functional equation 

a' 

An important symmetry of string models is the Target duality R' = — . Is there 

R 

an analogy of this symmetry in the Green's function expressions (4.74) or (4.75)? 
From the relation of the Green's function to the Lerch zeta function 

1 00 —x 
2*G(x,t;0,i*)\R = ^ £ + l)" a e ^ 

71=0 x (4.83) 

1 R„~ R i27rn— — Iff.* R, 



= p(-rE( n +-r' e 2nR = r(- 

hs n — ' n hf n 



R x a' z -l v " ' a' " R x a' ~ rx 2nR ,S ' a ] 

n=0 
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and the corresponding functional relation (2.12) we have: 

x 2itR 

1 /r, in * — 1 — * n 



e mx/R = (2iT)- s T(s)e ms / 2 e a £ ^— e a 

1 (4-84) 



in 

+(2ir)- s T(s)e- ms / 2 e a 27riT £ 2 e a 



These relations can be derived using a theta representation for the Green's function 
by considering the integral/ dC £ e a e 27ri? C s / 2_1 . Now setting: 

n=0 

2vri? 

1 • , D/ „ 1 -in- 



E e ~ inx/R ' = E x — e a 

"(„ + £L^ " (n + 



a 



2vri? y 

2vri? 

I in- 



inx/fl" = ^ e fl 



We obtain: 



For the case relevant to the zeta function x = ira and 

x' = 7ra 

x" = 2ttR - vra 
R" = a - 



a 2 1 



2 i? 



(4.85) 



t' = -t 

a' = a" = a — a 
x' = x 

, = ax 1_ (4.86) 

~ 27T.R 

x" = 2irR - x 
2ttR 



(4.87) 



If one sets R = a as in the previous section we have x' = x" = x = ira and 

R' = R" = R/2 = a/2 . If instead one sets R = —a then R = R' = , 

\/2 v 2 

x = x' = ira , x" = (y/2 — l)na and R" = (1 i=)a . Thus we find that R = -^= 

v 7 ^ V2 
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is a self dual point with respect to the transformation 

# = £± (4.88) 

In string theory the duality relation says that string propagation amplitudes on a 

a' 

circle of radius R can be written in terms of string propagation on R' = — . For the 

R 

Riemann point particle the duality reflected in the functional equation (2.12) tells 
us that point particle fermionic propagator from (x = 0, t = 0) to (x, t — ia = —isa) 
on a circle of radius R can be written as a superposition of the propagator from 
(x = 0, t = 0) to (V = x,t' — ia' = — s)a) on a circle of radius 

< 4 - 89 > 

and the propagator from (x = 0, i = 0) to (x" = 2-7ri? — x, t" — ia" = — s)a) on 
a circle of radius 

= 4 (4.90) 

Thus the radii, position and a values relevant to the Riemann zeta function R = 
a, x = ira, a = — are special values with respect to the duality transformations. 

4.6 a = i and path integral zeros 

Consider the integral representation of the zeta function where the integrand can 
be represented in terms of theta functions: 

S oo 

^ = 2 i-!_ 12 /_ 1 ^-^7 / ^t*/ 2 (0 4 (O|*t) + 9 2 (0\ir) - 9 3 (0\ir)) (4.91) 

If one can equate this integral with the result of a path integral of a quantum 
mechanical system then any statement about Riemann zeros can be translated to a 
statement about path integral zeros. 

Path integrals over gauge fields and fermions can develop zeros in certain situa- 
tions [21]. An especially simple example is + 1 dimensional gauge theory coupled 
to a single fermion with path integral [22], [23]: 

i 

DA D4;e° (4.92) 

The action has the gauge invariance ip'( T ) = e l "( T )^(r), A' q (t) = Aq(t) + 8q6{t). 
Note that this Aq(t) should not be confused with A\ parameter of section 3 as Aq(t) 
is a worldline gauge field and A\ is a target space parameter like R . Similarly in 
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this section ip(r) is a worldline fermion and ip(x, t) is a target space field. Fixing 
the gauge Aq (t) = ao for the + 1 dimensional gauge field then a global symmetry 
of the gauge fixed action is ao — > «o + 2ir. The path integral after gauge fixing and 
integration over the fermion reduces to: 



oo 

J da cos(a /2) (4.93) 



The path integral is not invariant ao — > a® + 2ir . There is a global anomaly. 
Indicative of this is the fact that the path integral vanishes. What's more is that 
the path integral with the insertion of a gauge invariant function vanishes[21]. That 
is : 

oo 

Z= J da cos(a /2)f(a ) = (4.94) 

— oo 

if /(ao + 27r) = /(ao) • However the path integral is invariant under ao — > — ao 
which is world line parity. There is a tension between parity and gauge invariance 
in this model. The gauge invariant form of the path integral is : 

Z> = J da -(e iao + 1) 

2w ° , (4.95) 

= / da e ia °/ 2 cos(ao/2) = / DA Q D^eh <MV>*(V>+^)+^ /2 



This gauge invariant expression Z' is invariant under ao — > ao + 2tt but not worldline 
parity. The parity invariant expression Z is invariant under ao — ► — ao but not gauge 
invariance. 

In our case we have a point particle Green's function which can be defined by 
the path integral: 

x(l)=x,t{l)=t 

G 4 r(x,t;0,0)\ R!A > = J DpDEDxDtDeoD^e- 1 (4.96) 

x(0)=0,t(0)=0 

Fixing the world line reprametrization symmetry by eo(r) = C we can use the 
methods of section 8.4 to express the path integral as: 

n + A\ 9 n + A'. 

00 ,c i 00 -C( ) if )x 

G*(x,t;0,0)\R,A> = I dCC*/ 2 - 1 E e R e R 

°° n + A' n + A' (4.97) 

n=— oo tt 



47 



Forming the superposition: 



Z =G 4 ,(x = 0,t-i(j;0,0)\ RtA/=0 

- G 4 '{x = 0,t- ia; 0, 0) | 1 

R,A'= g 

- G 4 , (x = nR,t- ia; 0, 0)|^ =0 (4.98) 

where we have introduced the a parameter through t — > t — ia as in section 4. 
In forming this combination the second term in (4.97) drops out and we have the 
simplified expression: 

Z = (2 s - 1)(2 X - S - 1)tt- s / 2 T(s/2)C(s) 

oo 

= J dC c s / 2 -He 3 (o\^)-e 2 (o\^)-e 4 (o\^)) (4.99) 


Introducing a Liouville type variable through C = a?e ip we have: 
00 222 

z = J d^){em a -^)-e2^\ a -^)-em a -^)) (4.100) 


To see why we call this a Liouville type variable consider the first term in (4.100) 
which is of the form: 

oo ^ a 2 n 2 

Zi = ( d<p(cos(tp/2) +ism(t( P /2))e ipa/2 ^ e" 6 # 2 

-oo n=-oo 

= /^ D P D ^/* lP±+< "* + r- e, ' ( " 2p2)+ '''* (it4 ' ) '" (4.101) 

In writing this expression as a path integral we have formed a superposition of 
antiperiodic (ip(0) = —tp(l)) and periodic (-0(0) = fermionic path integrals 

accounting for the cosine and sine terms above. The quantization of P follows from 
the periodicity of the x field. The Liouville interpretation follows from the fact that 
P 2 acts like a cosmological constant in + 1 dimensions. 

For the special case a = -a and using the modular transformation properties of 

^ 2 2 2 

the theta functions e^^Ole*^ )- ^(Ole^-^) - 64(0^^)) is an even function 
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of ip. The integral with the sine term in (4.101) is then identically zero and we have: 

Z = J ^(cos(^/2) e 4> 3 (0K-2) - 9 2 (0\e^) - e 4 (0\e*—)) = 

-oo R A R z R A (4.102) 

rn n P n n, So *-(p±+W+ \f ~ e*{a 2 P 2 ) + ^(itopW 
= J DipDPDxDipe 4 

At a zero of the zeta function so = to + i— we have a path integral similar to the 
example of the beginning of this section and we have: 

= Z(s ) = (2 s - 1X2 1 - - l)n~ s / 2 T(s/2)as)\ s=Uo/a+1/2 

= / ^(cos(^/2)e^/ 4 (0 3 (O|^) - WWjp) - OM^jp)) 

oo ^2 ~2 ~2 

= / da ( C os(a /2)e a o/ it o(e 3 (0\ie a o/ t o—) - 2 (o|i e <V*o ) - 6 4 (0\ie a °/ to -^)) 
= / DA DPDxD^eSo dr(Px+ri>+Ao/ito-e A o/to(a 2 P 2 )+r(iAoW 

(4.103) 

Where we have defined ao = io¥>- 




Figure 17: Integrand for the theta representation of the zeta function for the first 
Rieman zero 

In figure 17 and 18 we plot the real part of the integrand in (4.91) and its product 
with cos(toy/2) for the first nontrivial zero of the zeta function. Note like the gauge 
example at the beginning of this section it is the cosine factor which is crucial in 
obtaining a zero, the theta functions occur whenever the x field is periodic. In terms 
of the additional modular symmetry C — > — , the variables C = e v transform as 
if — > —if. This is worldline parity in the transformed variables A$ — > — Aq . In 
terms of the original wordline variables C = eo the usual modular transformation 
C — > — C yields worldline parity eo — ► — eo of the point particle einbein. Thus the 
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Figure 18: Integrand for the theta representation of the zeta function without the 
cosine factor. 

relation of worldline parity to the usual and additional modular transformations is 
reversed upon the introduction of the Liouville type variable <p. 
One can also use the integral representation 

Z = (1 - 2 1 - s )I»C(s) = J dCC s - x ^- x (4.104) 



to obtain a path integral representation of the zeta zero condition. In that case: 

00 1 
= Z(s ) = (1 - 2 1 - s )I»C(s) = / dCC 8 - 1 



00 00 



«T e^ + l 



= £ / dCC^- l C l / 2 e~ c ^ n+l ^ 

n=0 

00 00 _ C(n-\ — ) — nr(n H — ) 

= iJ2JdCC ° ~ x C x l 2 e 2 y 2> y 2 1 

n=0 

e ^ 2 W /2/ X N ■ / !x 

00 00 _ e *v (n H — ) — nr(n H — ) 

= z£ / ^cos(t (^/2)e^/ 4 e 2 y 2 J y 2> 

n=0— 00 

w e 1 ^/ 2 1 

. , n n n n . Jo (PQ+4>'i>+ T 9 (iff + e^ /2 )-(p 2 + <Z 2 ) - ito^V 
= 1 J DipDpDqDipe 4 2 z 

e A /2t Q ]_ 

/^^Mo/^-i^W-— = (i7r + e A «/ 2 'o)-(p 2 + (? 2 ) 

= z J DA$DpDqDipe 11 

(4.105) 

Here we have defined the Liouville type field by C = e^l 2 = e A °/ 2t °. The quanti- 
zation comes from the harmonic oscillator Hamiltonian -(p 2 + q 2 ) with eigenvalues 

(n + -) instead of from the periodicity condition. 
4.7 a as an extra dimension 

The transformation s' = 1 — s or cr' = a — a could be understood as a reflection or 
parity symmetry if a could be considered as an extra dimension. Also the spectral 
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properties of the Dirac operator can be investigated from the point of view of extra 
dimensions especially with respect to its chiral properties. Examples of this include 
fermionic field theory in the presence of a defect, overlap and domain wall approaches 
to lattice fermions [24], [25], [26]. 

To begin the investigation of a as an extra coordinate we write the Dirac equation 
equation (3.54) as: 

(a% + a 1 d x Wt,x) = {<J° — (e iadt -l) + a l d x )^{t,x) =0 (4.106) 

ai 



Where we have defined dt = —(e iadt — 1). Writing this equation in matrix form we 

la 

have: 

( ° _ D )(*) = ( _° * + M(lM=0 (4.107) 

\ —D J \ 1p2 J \dt-d x J \ ^ J 

Now we can modify the above equation by replacing the zero diagonal components 
by operators M and M so that the equation becomes 

I a D\f, Pl \(M B, + D.\/*) =0 (4.108) 
V -D M A *2 / \ »t-d x M ) \ i' 2 j 

The usual choice for M = ml leads to the massive Dirac equation whereas M = 
d a — msgn(cr) leads to the overlap equation with a mass acting nonuniformly in flavor 
space [26]. Neither of these choices leads to the zeta function as a Green's function. 

- - 1 _ 8 

Instead we choose M = d a + idt and M = d a — idt where dt = — (e * — 1) and 

B a = —( e ~ aid - - 1). Then we have: 
— ai 

( M . _^M(M =0 (4.109) 

y —D M J \ ip2 J \ d t -8 x d a -id t ) \ij 2 ) 

After Fourier transforming the above equation using the mode decomposition 
e -iEt+ipx+tku above equation becomes 



(e ak - 1) - (e 1 ^ - 1) {e ah -l)+pa \j Vi 

V (e aE -I) -pa -{e ak - 1) - (e iaE - 1) J \ V 2 

Searching for right moving solutions as above we set 



(4.110) 



(4.111) 
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and obtain the equations: 

(-(e ak -1) --(e iaE -1))^ =0 

f a n (4.H2) 

(-(e^-l)--)Vi=0 

Ob it 

Solving this equation we have E = - log(a— + 1) and k = iE so that the mode 

a K 

solutions are of the form 

x 

Y> n (x, t, a) = (^)- {it+a) (n + -)- lt ' a {n + -)-*/* e m R (4.113) 
and the Green's function is given by: 

G(x, t, a; x',t', a') = £ M*, *> O 

n 

. (x ~ x') (4.114) 

which again is proportional to the Lerch zeta function. 

In deriving the mode expansion one can also start with the 4D Dirac equation 
for massless fermion: 

(a°d t + a x d x + ct 2 <9 ct + a 3 d z )^(t, x, a, z) = (4.115) 

where we have used the Weyl representation 

\ 



a» 
J 



(4.116) 



In component form the corresponding equation is given by: 

/ a D\f^) = (s.-0. s l + a, w*\ o (4 . 117) 

V -D M ) \ f. 2 ) \ 8,-8, a, + id, J\tp2 J 
If one identifies z coordinate with t we obtain an equation in the form (4.110). 

4.8 a and Space-time Reflections 

The transformation s' = 1 — s, s" = s* ,s"' = (1 — s)* of the Riemann zeta function 
have interpretations as discrete space-time reflections when a is interpreted as an 
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extra dimension. Using s = 



it + a 



we have 



a 



it' + a' = it + (a-a) 

it" + a" = -it + a 

it'" + a'" = -it + (a-a) 



(4.118) 



and thus the time reversal T and parity P transformations are given by : 



T : t 



-t 



a — > a 



PT : t — > —t, a — > a — a 
P :t t, a — > a — a 



(4.119) 



Note that the parity transformation includes a small translation by a in its definition. 
This reduces to the usual parity transformation at large distances a » a. 

The transformation properties of the Riemann zeta function under s — > s* and 
s — > 1 — s have well known applications to the Riemann zeros. If s is a zero than 
so is s* ,1 — s and 1 — s* . This means that if a zero in the critical region exists 
there would actually be four zeros. Thus from the transformations (4.120) for P,T 

and PT one can restrict the search for zeros to the region — < — < 1 and t > 0. In 

2 a 

the interpretation of a as an extra dimension the zeros have the interpretation as 
locations (x, t, a) that thatfor which the fermion has zero probability amplitude to 
evolve to, essentially forbidden points of transition. 

5 Statistical mechanics of the fermionic theory 
5.1 Partition function 

It is known that the zeta function behavior at s = 1 can be interpreted as a thermal 
divergence [27], [28], [29]. Also the mixing parameter a parameterizes a Boltzmann 
distributed initial state so it is interesting to study the Weyl fermion theory at 
finite temperature to see how this phenomena is represented. In this section we 
compute the finite temperature partition function associated with cases (2) and (4) 
with Hamiltonians 




(5.1) 



and 




(5.2) 



respectively. The partition function is defined by: 



Z(f5) = tr(e-^ 2nd ) 



(5.3) 
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with i?2nd the second quantized Hamiltonian. Because we use the second quantized 
Hamiltonian we will compute the grand canonical partition function which contains 
arbitrary numbers of particles. Nevertheless we shall show that it is the one parti- 
cle partition function z one (P) which controls the behavior of the partition function 
near the thermal divergence. The general formula we shall need is that for a non- 
interacting fermionic field theory the logarithm of the partition function is given 
by: 

oo 

log Z = !og(l + e-ft en -ri) (5.4) 

n=0 

where j3 is the inverse temperature and fx is the chemical potential. 
Case (2) e n = 2-nn/L 

Again we begin by reviewing the well studied case of a free right moving fermion on 
a cylinder of circumference L = 2irR with Hamiltonian 



jdvitW (5.5) 



H 

In terms of Fourier modes it can be written: 

oo 

H = ^ e n b- n b n = — XI nb^ n b n (5.6) 

So that: 



2vr - 

n=l n=l 



Z(P) = f[(l + e 27rn L) (5.7) 

71=1 

If we had studied a scalar particle the partition function would be: 



_ 2nn— , 

^scalar(/3) = II (! " e L ) ( 5 ' 8 ) 
n=l 



In appendix A we discuss the treatment of the scalar particle associated with the 
nonstandard dispersion relation (3.31). 

For the fermionic theory the number of states at energy E is given by the inverse 
Laplace transform: 

ioo 

a(E)= f dpe^ E Z(P) (5.9) 



-j^ oo 

As E = — X nN n then a(E) gives the number of ways of writing £ as a sum of 



R , 

71=1 



positive integers irrespective of order with no repeats. It is this property of the 
partition function that plays a role in additive number theory where one seeks to 
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enumerate the ways of writing a positive integer as a sum. 



Case (4) ae n = log(l + a2im/L) 

In this case we have H = J dxtp — \og{iad x + l)ip with periodic boundary conditions 

x 

oo —ie n t+in— 

ip(x + 2ttR) = ip(x). Fourier transforming ip{x) = J2 b n e and using the 

n=0 

j ji 

dispersion relation for the one particle energies E n — - log(a— + 1) we have: 

a R 

°° 1 271-r? 

H = V-log(a— n + l)6_ ra 6„ (5.10) 

CL Li 

n=l 

Then from the formula for Fermi-Dirac partition functions we obtain: 

oo oo 

z(d) = n a + = n (i + + v- p,a ) ( 5 - n ) 

n=l n=l 

A general result on infinite products states that an infinite product such as 

oo oo 

n (1 + o-n) diverges whenever J2 a n diverges [30]. Applying this to our case we 

n=0 n=0 

have 

oo 

n (l + e _ ^ £n ) diverges whenever the one particle partition function z onc ((3) = 

n=0 

oo 

e -P £ ™ diverges. Thus the grand partition function diverges whenever the one 

n=0 

particle partition function diverges. For the dispersion relation (3.31) the one par- 
ticle partition function is 

*™o?) = {-f /a £ — ^— = (-r /a tH^, -) (5.i2) 

° ^o(n + -)^ a a 

which has a simple pole divergence only at (3 = a. Thus we see that the grand par- 
tition function of the nonstandard fermion also has a divergence at (3 = a. Studying 
the Yang-Lee zeros of the one particle partition function yields another physical 
representation of the Riemann hypothesis. We discuss the one particle partition 
function in Appendix B. 

Now consider the fermionic partition function associated with the special case 

oo 

R = a . The formula for the total energy is then E tot3i \ = - ^ (l°g n )An and the 
partition function becomes: 



a 

n=2 



Z(J3, L = 2ita) = H (1 + n-< 3 l a ) (5.13) 

n=2 

From the definition of the partition function the number of states at a given energy 
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is given by: 



a(E) 



dPe 0E Z(l3) 



(5.14) 



The usual situation is that the number of states grows with the energy. However 
because of the definition of energy a mu \ t (E) gives the number of factors of E without 
regards to order and with no repeats. For example E = log(120)a _1 the states 
without regards to order with no repeats are of the form FJ &-n |0) 

N n =0,l;n=2,...,exp(aE) 

and can be listed as: 

120 = 5! = 

5- 4-3-2;&_5&_4&- 3 M 0) 
20-3-2; &_ 20 &-3&- 2 | 0) 
15-4-2;6_ 15 6_4&- 2 | 0) 
12-5-2;6_ 12 6_ 5 6_ 2 | 0) 

10 ■ 4- 3;6_ 10 6_4&-3| 0) 
8-5-3;&_ 8 &_ 5 6_ 3 | 0) 

6- 5-4;6_ 6 6_ 5 6_ 4 | 0) 

(5.15) 



60 


2; &„ 60 6- 


-2 


0) 


40 


3; 6_ 40 fc- 


-3 


0) 


30 


4; b- 30 b. 


-4 


0) 


24 


5; 6_24&- 


-5 


0) 


20 


6; &_2o&- 


-6 


0) 


15 


8;6_ 15 6- 


-8 


0) 



12 



10;fe„i 2 6 



120-10 



0) 



120;6_ 120 | 0) 



There a(- log(120)) = 15 states without regard to order with no repeats. These 
consist of 1 four-particle state, 6 three-particle states 7 two-particle states and 1 
one-particle state. A dramatic effect occurs when E is the logarithm of a prime 
number in units of a. Then there is only one state and <J mu it{E) = 1. This behavior 
leads to an inhomogeneous density of states and a potential loss of equilibrium. 

Returning to the general case of arbitrary R = L/2tt the partition function can 
be written: 



27T71 



Z(J3,L)=H(l + (a— + l) 



(5.16) 



For L » a we obtain the additive partition function (5.7). For L = 2ira we obtain 
the multiplicative partition function with a mu it(E = log(Prime)) = 1 . For the 
additive density of states the number of states grows rapidly with energy <r a dd(-E) ~ 
e 2 V wLE / 24 = e n V RE / 3 . The general expression a(E, L) interpolates between these 
two extremes. 
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5.2 Series representation of the free energy 

From the representation (5.4) specialized to e n = - log (a— + 1) we have: 

a Ix 

oo 

logZ = ]T log(l + + l)-^°e^) (5.17) 

oo (-y) m 

Now expanding the logarithm in the sum using log(l + y) = — J2 we obtain: 

m=i rn 



(5.18) 



logZ= E E -(-l) m -(a- + l)-^ m / a e^ m ) 

m=ln=0 mi? 
= ^ ( l) m 1 ( R \ m l 3 / a c l 3 ^ m ^ H (^ m ; ^) 

m= i ma a ' a 

For the special case of = o and /x = this becomes: 

00 1 R 

logZ(R = a, f i = 0)=J2-(-l) m -a^ L ) (5-19) 

m=l 

In this form we see that the divergence at (3 = a is contained in the first term of the 
series. 

5.3 Thermodynamics quantities 

Various thermodynamic quantities can be obtained from the partition function. For 
example the free energy / , average energy U , average pressure P and average 
occupation number N are given by: 

1 1 - n ± x ± 

f((3,L) = logZ(/3) = -5>g(l + (a- + l) a)~-£n a 

' ' n=l ' ' 



n 



d 1 00^ log(a- + 1) x 



8(3 a^H P 



n=i n 



l + (o- + l)a 



d 1 d 

dV U = 2^dR U 



1 1 00 



1 + (0-5 + 1)0 



1 1 °° 



I 

!og(«^ + 1) 

Tl — — 

yl + (a- + l) a 



- 1 
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n = -£-EMi + (a; + ir /3/a e- a ) = E 



n=0 n=U ^ — 

e a + (a- + l)a 
H 

~ E n ^ /a ( 5 - 2 °) 
All these quantities diverge at the critical inverse temperature j3 c = a. 

5.4 Thermal Green's function 

The thermal Green's function for a fermionic field theory is defined by: 

CO 

Gp(x,t;0,0) =tr{e- pH ^{x,t)i){0)) = ^ (-^) m G(x, t + if3m; 0, 0) (5.21) 

m=— oo 

where G(x, t; 0, 0) is the zero temperature fermionic Green's function. From the 
definition we see that the thermal Green's function is antiperiodic in imaginary 
time with period 27T/3 . Now using the formula for the zero temperature Green's 
function associated with case (4): 

~, n 
1 i— X 

G(M;0,0) = J2 n t/a eR 
n=o (a^ + 1J 

,R^,„^ l ^ 



(-) it/a V e R' 

a 

,R^+/ n , , x it R s 



-o (n + -f 



= (-^X^,-,-) (5-22) 
a 27rit a a 

the thermal Green's function is given by: 

G /J (x,t;0,0)= f; (-ir(^)-^/>(^,^-^,^) (5.23) 

m=— oo 

For the special case relevant to the Riemann zeta function we have R = a, x = 0, t = 
so that 

oo a 

G^(0,0;0,0)= E (-l) m C(^) (5.24) 

m=— oo 

As the Riemann zeta function has a pole at z = 1 we see that the terms of the 

sum for the thermal Green's function develop a divergence whenever (3 = — and 

m 

m 

temperature kT c = — . 

a 

The crucial difference between the thermal Green's function and the mixing 
Green's function used to describe the Riemann hypothesis is that a trace is taken 
in (5.21). It is simple matter to introduce the mixing parameter a into the thermal 
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Green's function and one has: 

Gp( X ,t;0,ia)= £ (-ir(-Y U - mf3)/a <f>(^ + " ^ , ~) (5-25) 

m=—oc 

Again specializing to the case relevant tot the Riemann zeta function R = a,x = 
0, t = we obtain: 

OO n 

G (3 (0,0;0,ia)= £ (-l) m C(^ + -) (5.26) 

a a 

m=— oo 

In the presence of the mixing parameter a the terms of the sum develop a thermal 

divergence at (3 = - — — thus the temperature of the thermal divergence is modified 
m 

m 

to kT c = . 



a — a 



6 Harmonic oscillator representation of the zeta func- 
tion, 2D fermionic string theory and Matrix models. 

6.1 Logarithmic oscillator and fermionic string theory 

In the previous sections we obtained the zeta function as a Green's function from a 
fermion theory compactified on a circle R = a and described by the actions: 

1=1 dtdx{ii)*d t 4> + 4>*-log(iad x + 1)V>) (6.1) 

J CI 

or 

1 = J dtdx(i)*\e aidt - l)dt*l> + 4>*idxi>) (6.2) 
The sum over integers in the mode sum for zeta arose because of the compactified 

Tl 

momentum condition p = — and the nonstandard energy momentum dispersion 

11 

relation E = - log(ap + 1) or p = -(e aE - 1). 
a a 
Another method of introducing the integer n into the mode expansion is to use 

the quantization of a nonstandard oscillator with energy 

E = -\og(a\{p 2 +u 2 \ 2 ) + l) (6.3) 

CI — 

Then the Schrodinger equation for this system is: 

idtfp + - log(4(-df + lo 2 \ 2 ) + 1)V = (6.4) 

Ob Z. 

with energy quantization condition E = — log(ao;(n+-) + l) = — \og(aujn+-auj+l). 

a z a z 

The relation with the Lerch zeta function follows immediately upon forming the one 
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particle partition function 

1 1 

-8— log ( awn H — auj + 1) — na 

^io g -osc(/3,M) = E e-^"^ = E e a 2 

n=0 n=0 

(6.5) 

= H a E r^-j e~» n = (auj)~a ( p(iuAi — + \) 

aw 2 

Here we have introduced a chemical potential associated with oscillator number, 
^log -osc(/5, /i) is proportional to the Riemann zeta function for auj = 2 and u = 
and the Dirchlet eta function for auj = 2 and = m. The condition aw = 2 takes 
the place of the condition R = a of the previous section. The Riemann hypothesis 
can be formulated in by saying that the one particle partition function only has 
zeros for Im/3 = — when fx = m or that there are no phase transitions in the form 

of Yang-Lee zeros in the nonstandard quantum oscillator for -a < Im/3 < -. 

2 a 

The dynamics of the classical logarithmic nonstandard oscillator are given by 
the equations where (q = A): 



. , 8E V 1 P r,P 

Velocity : q = — = ^ § = ^e~ cL 

op m , v 1 9 9n m 



dE 

— = -mu"q 2 

ag p z 1 

2m 2 y ; (6.6) 



Force : p = — -777- = —muP'q ^ = —m^ 2 Qe aE 

a( 1 — mw g J + 1 



Position: g(t) = go cos(u>te ) -\ sin(wte 

2vr mUJ 
Period : T = — e aE 

W 2 2 

^ „*? / Pn mw 9 , 

Dispersion : e ah = a{-^- + — — gg) + 1 

For initial momentum po = we see that the period 

T=*U- = *V=£,8 + 1) (6.7) 

So unlike the standard oscillator the period is amplitude dependent. 

Returning to the Schrodinger equation (6.4) we see that this equation follows 
from the two dimensional field theory: 



I 



= I dtd\(iij*d t ip + V*- log(<4(-8j| + w 2 A 2 ) + l)ip (6.8) 
J a £ 



In the limit of aui — > this reduces to the fermionic formulation of the c=l matrix 
model [31]: 
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/ = / dtd\W*dtil> + r\{-d\ + ^A 2 )V>) (6.9) 
J la' 



1 ,2 



for — = u . In the c=l matrix model one can go from the right side up to inverted 
a' 

harmonic oscillator potential by sending a' — ► —a' [32]. When one does so one 
changes the background for the model from as well as the vacuum energy. One also 
replaces the Hermite polynomials by parabolic cylinder functions and the physical 
description is quite different in the two cases. 

The Green's function associated to the action is given by: 

G(t,X;t',X') = J dpe ipt J dse~ sp+i ^ s ^Ai|e iah |A 2 ) (6.10) 

with 

(\i\e iTH |A 2 ) =Y,e- iEnT H n (\)H n (\') = £ (ow(n + \) + l)- lT ' a H n (X)H n (X') 

n n 

(6.11) 

Unlike the case (4) of the previous section the Green's function does not yield the 
zeta function. However the one particle partition function defined by z osc ((3) = 
J dXG(X, X, (3) does yield the zeta function as we saw above. 

One of the many fascinating aspects of the usual c=l fermionic theory is that 
it is quadratic in fields yet reproduces the genus expansion for 2D interacting c=l 
noncritical string theory through the formula [34] : 

oo oo 

<^=\§j^j^<A A) = i^j^r J d\G(\, A; r) = 

-oo 

oo oo 

= / > = >4 Re *4 + #> 

(6.12) 

The genus expansion results from expanding the above expression in powers of /i 2ff_2 
where g is the genus. 

Here we study the modification to the above formula resulting from using the 
logarithmic fermionic oscillator theory (6.8) . From the structure of (6.12) we see 
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that we have: 



oo oo 

d^p = — Im J dre'^T J d\G(X, A; r) = 

-oo 
oo oo 

= ^-Im / dre~^Tz log ^ osc (a,r) = ±-lm f dre^r^r^ (-, — + h 
Z7T J Z7T J a au> 2 




i \ ) 

n -log(a(n + ~)0J + 1) +ifl 
a £ 



(6.13) 



Also for the special value auo = 2 the integrand becomes: 

oo 

VW = Eflm / dre-^r(aa;)-/ a C(-) 



1 

= E T^Im / dre-^ T T{aLor T/a C(.-) (6-14) 



oo 

+ E^Im [ dTe-^ T T{aio)- T / a TT (l-p-^)" 1 

m 27T J 

peprimes 

where the first integral is over the critical strip and second integrand involves the 
Euler product over prime numbers. It would be interesting to investigate if this 
occurrence of prime numbers has a connection with the adelic 2D string theory 
studied in [35]. 

Another manifestation of the c = 1 2D string model is the compactification of 
the t direction with radius R . In the usual case we have [36]: 

2ttR 

1= J dt J dX(iij*d t ij + ij*^(-dl+Lo 2 X 2 ^) (6.15) 
o 

with the result: 

d llP (p,u;,R)=d^Re(Y / j ^ ) (6.16) 

mn (n + -)oj + 1 + (m + -)— + ifx 
2 2 ti 

We compute the modifications the take place with the logarithmic oscillator action 
(6.8) 

2-kR 

1= J dt J dX{ii(j*dt^ + i)*^\og{a^{-di+Lo 2 X 2 ) + l)ilj) (6.17) 
o 

As in the usual case the propagator is antiperiodic 
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G r {X,X';t) = £(-l) n G(A,A';r + 2^Rn) an d we have 

n 

oo oo 

8^p= i-Im y dre^r | cIXGr(X, A;t) 

-co 

oo 

= E T^Im / dre-^ r rz log _ osc (a,r + 2vrmi?)(-l) 

m 27T J 





oo 



(6.18) 

= E Tr Im / dTe ' lflT ■ u"/9^ r W' T/ °Cg(- — + J) 

m 27T J smh(r/2it) a aw 2 

o 

= ^ Re (£ i 1 n ) 

™« - log(a(ra + -)u + 1) + (m + -)— + i/x 
a z z ri 

In computing the above we have used the fact that the finite R dependence can be 

implemented by inserting the factor f— into the integrand as in [36]. This 

sinh(T/2i?) 

expression agrees with the usual vacuum energy (6.16) in the limit a — > . Note 

a' 1 

that the duality symmetry R — > — = — ^ — of (6.16) is violated by the a dependence 

R u> R 

in this Harmonic oscillator representation. Previously we found a generalization of 
the duality symmetry using the Green's function representation in section 4. 

6.2 Matrix Model 

It is well known that the fermionic action (6.9) has an interpretation as a matrix 
model [36]. Does the nonstandard fermionic action (6.8) also lead to a random 
Matrix interpretation? If we write the matrix M(t) in terms of it is eigenvalues 



/ Ai(t) ... \ 
M(t)= : ■.. : 

V • • • X N (t) J 



(6.19) 



The fermionic action is associated with the Matrix integral 

, / dt(tr(M(t)n(t)- — Ltrlog fv / a 7 (n(t) 2 + ^M(t)) 2 + 1) . 

Z= DMDUe Va' \ a J ( 6 2 o) 

where II(i) is a conjugate auxiliary field. In the eigenvalue basis the integral be- 
comes: 

f fdt(\p logfCt hi) 

Z= DXDpA(X)e a SV 2 ' (6.21) 

with Vandermonde determinant A(A) = Jl (A(i) — X(t')) arises from the measure 

of hermitean matrices . It is the antisymmetry of the Vandermond determinant 
which connects the matrix integral with the fermionic theory [31]. 
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The eigenvalue action is somewhat more complicated than is usually considered 
for Matrix models. After eliminating the auxiliary field p we find the eigenvalue 
action: 

S = J dt(Q(X, A) - V(A)) (6.22) 

with 

Q(A, A) = A 2 2 (1 + ^f) 

l + ^/l-2aA 2 (l + ^) 

--log( ■ 2 ) (6.23) 

l + ^l-2aA'(l + ^) 

nA) = ilog(l + ^) 

A 2 A 2 

This reduces to the usual matrix eigenvalue action f dt( uj 2 — ) in the limit of 

auj — > . In figure 20 we plot the potential F for the Matrix model. One can add 
an inverted oscillator potential to this matrix action: 

V(X) = -^X 2 + - a log(l + a^A 2 ) (6.24) 

This form of the potential has been shown to survive large N scaling and has been 
used to describe a Dq brane contribution to type OB string models in two dimensions 
[37] . Shown in figure 19 and 20 is the potential energy used in that model. Taken 
with the analysis above Dirichlet series may play a fundamental role in the short 
distance dynamics of string models. 




1 auj 2 X 2 

Figure 19: Logarithmic harmonic oscillator potential V(X) = -log(l H — ). 

a — 

Setting (3 = Ne one can write the path integral representation for the ordinary 
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V(A) 




Figure 20: Sum of inverted and logarithmic potential. The same potential was used 
to describe OB string theory in two dimensions in [37] 



oscillator as 



r, (a \ rn a a + ~T q ^ ~ ipnQn + ^n+l , 

Z \P) = J II • • • dq n dp n e 2 dg n _i . . . 

n=l 

oo 

= / II • • • dq n dp n (p n -i\ q n ) (q n \T\ p n ) (p n \ q n+ i) ■ ■ ■ 

n=l 

oo 

= / II ■■■dq n (q n -i\T\ q n ) (q n \T \ q n +i) ■■■ = Tr{T N ) 

n=l 



(6.25) 



where the transfer operator in the mixed basis is given by: 

I 2 , 2 2\ 

{q \T\ p) = e- iPq - £ 2 iP +UJ<1) 



(6.26) 



The partition function for the ordinary oscillator is: 

^(/3) = E^ (n+1/2)u = e "^ T 1 

71=0 



e-P" 2sinh(/3u;/2) 



(6.27) 



For small (3 this is z osc ((3) = ——. 

For the logarithmic oscillator one can proceed in a similar way as: 

U> 2 2 

„. a . , °° , , -£log(a( p 2+ — q n ) + 1) - ip n q n + ip n q n+1 
Z (P) = J 11 ...dq n dp n e 2 dq n ^i . . . 

n=l 

oo 

= J U. ■■■ dq n dp n {Pn-l\ Qn) (qn\T\ p n ) (j>n\ q n +l) ■ ■ ■ 
n=l 

oo 

= / 11 • • • dq n (q n -i\T\ q n ) {q n \T \ q n+1 ) ...= Tr(T ) 
n=l 

(6.28) 

So that in the mixed basis the transfer operator is: 



e-log(a^(p 2 +iu 2 q 2 ) + l 



(q\T\ p)=e a 



e -iW( a I(p2 +w 2 g 2) + 1 )- e /o (g^g) 
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The partition function for the logarithmic oscillator is then: 

00 B 1 1 

^io g -osc(/5) = E + V2)w + l)^ /a = M^Ch (S ? + — ) (6.30) 

n=0 

Near (3 C = a this is zi og _ sc(/3) = (aw)"^" 



/3-a 

6.3 Zeta function in the Topological Matrix Model 

Another matrix model related to 2d string theory is the Penner or topological matrix 
model used to compute the Euler characteristic of the moduli space of Riemann 
surfaces [39], [40]. The partition function in this case is given by the single matrix 
integral: 

Z(A) = J dMe -^ MA ^ v ^ tr = J dM(det M ) v - N e- vt < MA ^ (6.31) 

This matrix integral is a matrix analog of the Gamma function. As the contour rep- 
resentation of the Gamma function is so similar to the contour representation of the 
zeta function it is worthwhile investigating whether there is a Matrix representation 
of the zeta function as well. Indeed one can introduce the matrix integral: 



\v—N ^—mvtrMA 



Z\A) = [ dM(det(M)r- N eUt J A _ i = W dM(det(M)Y 

m=l 

(6.32) 

Now if we scale M by m~ l and take into account the scale transformation properties 
of the measure together with the Vandermonde determinant we have: 



- - - ■ - SN u—N „—mvtrMA 



Z'(A)= £ fdM(det(M))"-»e 

m_1 uN-NZ+NZ+ — { — 



(6.33) 



Thus we can express the zeta function as a Matrix integral through 

COO = z^jf dM(det(M)y- N eU J A _ i (6.34) 
N(N - 1) 



where s = vN + 
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6.4 — = - inverted oscillator and Riemann zeros 

a 2 

Recall the relation between the inverted oscillator and it's density of states [41]. 
The Hamiltonian of the inverted oscillator is: 

H = -{p 2 — uJ 2 q 2 ) = —-(x + X- + X-X + ) (6.35) 
where x± = (ujq ±p)/y/2 . A simple set of wave functions are given by [41]: 

oo 

iJj E (x + ) = ( dk(k/u)- iE l"- a ' a e ikx + = fe«y«+*/«-ir(l - - - %-) (6.36) 
J i a u> 

o 

These are energy eigenstates for — = - so that: 

a 2 

1 

x , iE/ui- - 1 E 

4> E (x + ) = fy 2F(--z-) (6.37) 
Placing a zero condition at q = A" 1 enforces a quantization condition: 

0(E, x+) = Imlog(^ B (s+)) = (n + -)tt (6.38) 

The density of states is given by p(E) = — = — (— — log A). 

dE Z7T dE 

Here we note the similarity between (6.36) and case 3 defined by the dispersion 

relation E = — log(p/a) with wave function: 
a 

oo 

Mx) = f dp{pa)- u / a - a / a e ipx = (-f/o+v/^Til - it/a + a/a) (6.39) 

J (X% 


The correspondence is: 

a/a -► 1/2 

t/a E/uj (6.40) 
x — > x + 

In section 8 we show that one can go from expressions like (6.39) to the zeta function 
by the process of periodization. To see how this works in this case form: 

Mx)\ R = £ C'^ f'^I^Til-tt/a + a/a) (6.41) 

m=— oo 
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Using the definition of the Hurwitz zeta function this can be expressed as: 



a 



+(-i)-iC ff (i-5,i-^))< 1 -r(i- a ) 

with s = (it + a)/a . If we choose x = nR and R = a and use (2 1_s — 1)C(1 — s ) = 
Cij(l — s, 1/2) we have: 

Mx = ™)\ R=a = (2n) s ' 1 C(l - s)(l + (-l) s - 1 )(2 1 ^ - l)^- s r(l - s) (6.43) 

Now using the functional equation: (27r) s_1 r(l — s)((l — s) = . — - this be- 

2 sm(s7r/2) 

comes: 

Mx = ™)\ R=a = -((s)(i s - rT 1 ^ + (-i) s ' 1 )(2 1 - s - i)r* 

= -cm-iy - 1)^(1 - (-i) s )(2 1 - s - 1) = c( s )(2 1 - s - 1) = -v(s) ' ! 

The formula that relates the phase of the zeta function at a /a = 1/2 to the 
number of zeros on the positive t axis to a value of T is given by [6]: 



2vr J R / (6.42) 



N(T) = -0(T) + 1 + -Im / (logC(s)) f ds 
it i Jc 



I , l c (6-45) 

= _^(T) + l + -Im(logCb+iT)) 

7T 7T 2 

it 1 t 

where $(t) = Im(logr(— + -)) — -ir and the contour C is the boundary of the 
rectangle defined by— e < Res < 1 + e, < Ims < T. The formula is not modified by 
the 2 1_s — 1 factor in (6.44) as both ((s) and r](s) have the same number of zeros 
in the critical strip. Formula (6.45) means that the process of periodization and the 
correspondence (6.40) to the quantum inverted oscillator leads to a condition on the 
number of zeros in place of the quantization condition (6.38). If we define: 

^ E {x + = tt) = r(i| + \)i~ E C{\ + iE) (6.46) 

Then the zero condition after periodization becomes: 

ttN(E) = Im(log(^(x+ = tt)) (6.47) 

Here N(E) is the number of zeros in the function /(e) = ip e (%+ = tt) up to energy 
E which corresponds with the Riemann zeros. It is somewhat surprising that pe- 
riodization through (6.41) leads to such a relation as the usual derivation of (6.38) 
uses the asymptotic form of the wave function at oo and periodicity does not permit 
such a region. 

Various other permutations of the dynamical variables can be considered besides 
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t/a — > E/oj and these are considered in appendix A. Often the physical picture 
changes including the introduction of a nontrivial background field due to the mod- 
ification of the equations of motion under the permutation. For example, if instead 
of x — ► x+ one chooses t/a — ► E'/o; and x — > p then the relation t = — {xp + px) 

becomes E = ^(xp + px). This is the Hamiltonian of Berry and Keating [2] as well 

as the dilatation operator. If one uses the original variables t = ^(xp + px) and 

H = -log(ap) then t has the intepretaion of the time of arrival operator. Usually 
one has difficulty in defining a time operator in quantum mechanics because time 
takes its values on the real line whereas energies are positive. It is for this reason that 
the momentum-position uncertainty relations and energy-time uncertaity relations 
have very different derivations. However in this case because of the the logarithmic 
dispersion relation the energy takes its values on the real line the treatment of time 
as a quantum operator is less problematical. 

7 Ground rings, algebraic curves, modular forms and 
Dirchlet L-series 

There is a remarkable connection between the inverted oscillator with t periodic 
with period (3 C = ^fa! = — and topological string theories on non compact Calabi- 



Yau manifolds [42] which are determined by algebraic equations. If this connection 
generalizes to the above process of periodization of x + with period R c = a it would 
represent a bridge between the zeta function, 2d gravity (quantum surfaces) and 
the algebraic geometry of complex manifolds. We investigate some of these issues 
in this section. 

7.1 Elliptic curves 

Projective elliptic curves E are defined as the points (X, Y) satisfying the equation 



The affine version of the elliptic curve is obtained by setting Z = 1. Elliptic curves 
are characterized by two quantities. The discriminant which is given by: 



[43]: 



Y 2 Z = X' 3 + aXZ 2 + bZ 3 



(7.1) 



A = -16(4a 3 + 27b 2 ) 



(7.2) 



and the j invariant which is defined by: 



3 = 



1728 • 4a 3 
A 



(7.3) 
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When there is a term on the right hand side of the elliptic equation of the form 
cX 2 Z these formulas are modified by (a — > a — c 2 /3) and (6^6 + — c 3 — -ac) . 
The discriminant may also be written as A = 16(ri — T2) 2 (r2 — r^) 2 (r\ — r^) 2 where 
r±,r2 and r% are the three roots of the right hand side of the cubic equation. 

The most famous elliptic curve doesn't actually exist [44]. Consider the curve: 

Y 2 Z = X 3 - (S + T)X 2 Z + (ST)XZ 2 (7.4) 

If s = r n and T = k n the discrimant of the curve is A = 16S 2 T 2 (S - T) 2 = 
16r 2n (k n — r n ) 2 k 2n . If k n — r n = i n for some integer i then there is no modular form 
associated with elliptic curve in contradiction with the Taniyama-Shimura theorem 
[45]. Thus i n + r n ^ k n which is Fermat's theorem. An algebraic equation of the 
form (7.4) is called a Frey curve. 

e~ aEk _ , 1 
Recall the Boltzmann probabilities Pk = — r- = {k a )t7 — considered in 

z{a) Cw<») 

section 4. If — is an integer n > 2 then besides obeying the usual identities < 1, 

a 

oo 

Pi + P r < 1, PjPj < 1 and P r = 1 the probabilities also obey the inequality 

r=l 

— / — + — because of Fermat i n + r n / k n . Rewritting this inequality as 



2P k {P l + P r )^2P i P r (7.5) 

This expression has the interpretation in terms of logical operations as P(kA(iVr)) / 
P(iAr). Thus the probability of choosing a pair of numbers a and b from a Boltzman 
distribution of the form above is never equal to the probability of choosing a number 
c together with (a or b) . To illustrate this consider the simple situation of three 
objects a, b and c with probabilities all equal to -. The 9 possible ways of choosing 
the two objects at random from {a, b, c} are 



^ aa ab ac ^ 

ba bb be 
y cc cb cc J 



(7.6) 



As we are assuming uniform probabilities in this simple example there are two ways 

2 11 

of choosing a and b namely (ab, ba) so the probability P(a A 6) is — = 2— • — . There 

9 3 3 

are four ways choosing c and (a or b ) namely 

(ac,ca,bc,cb) (7-7) 
4 111 

so the probability P(c A (a V b)) is - = 2- ■ (- + -) thus the inequality is obeyed in 

y o o o 

this case. 
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er aEk --log(a-| +A ) 1 ,-r, 1 
Returning to the probabilities P]~ = = e a R 



z z C( n ) 

we see from the discussion of section 3 that the fact that the k is an integer follows 

from quantum mechanics and the fact that wave function obey ip(x + 2irR) = ip(x) 

k 

on a periodic space, so that momentum is quantized p = — . The special nature of 
R R 

— = 1 and — = n follow from the duality symmetry (4.88) and the thermal partition 
a a 

function (5.11). For algebraic curves this means that if one considers an elliptic 
curve over a set of operators (like e~ aH considered above) which are quantized, 
than coefficients of the curve which would normally takes values in the complex 
numbers C or reals R become curves over the rationals Q. It is elliptic curves over 
Q which are in one to one correspondence with modular forms and Dirichlet L-series 
which in turn have a form of the Riemann hypothesis [43] associated to them. 

7.2 Ground rings 

In the Liouville theory of noncritical (c = 1) strings ground rings are defined from 
the operators [46]: 

x = 1 i(z) = (cb + idX-)e iX+ 

2' 2 ^ . y (7.8) 

y = 1 i (z) = (cb - idX+)e~ lX 

2'~2 

Where X ± = —j={X ± i<f>) an c and b are ghost fields. Defining the operators 

y 2 

X = xx,Y = yy, S = xy,T = yx yields four generators of the ground ring with one 
relation between them: 

XY - ST = (7.9) 

It is the relation (7.9) that corresponds to an algebraic surface, in this case the 
quadratic equation describing the conifold. 

In the matrix formulation of noncritical c = 1 string theory the ground ring is 
defined by the operators [46]: 



Oi )2 = (qu +p)e~ tu} 
2 ,i = {qoJ - p)e^ 



(7.10) 



With one relation between them 01,202,1 = ^ 2 q 2 — p 2 = —2H where H is the 
Hamiltonian of the inverted oscillator. 

For c = 1 noncritical string theory defined at A — D — E points at special radii 
R = n/(uj\^2) much more complicated algebraic equations are derived. For example 
for the Z?2 (n = 2) string theory the chiral ground ring operators where determined 
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in [47] and are given by: 

T=i(z 4 + y 4 ), Y= l -xy(x A -y% X = (xy) 2 (7.11) 

There is one relation between the ground ring operators which is: 

Y 2 Z = X 3 - (T 2 )XZ 2 (7.12) 

Here we interpret T as a Neron parameter associated to the elliptic curve E , Z is 
a parameter of the projective curve,. Again one can set Z = 1 to obtain the affine 
curve. 

The discriminant and j invariant of this curve are 

A = -16(4a 3 + 276 2 ) = -16(4(-T 6 ) = 64T 6 
= 1728 • 16a3 = i?28 (7.13) 
A 

As another example the non chiral D2 ground ring of c = 1 string theory at 
radius R = y/2a' which was also computed in [47] . The ground ring is given by the 
operators: 

X = xxyy 

Y = j(xx - yy)((xy) 2 - (yx) 2 ) 

r (7.14) 

T = -((xy) 2 + (yx) 2 ) 
{ 

S = -(xx + yy) 

Again there is one relation between these operators and it is given by: 

Y 2 Z = X 3 — X(T 2 + S 2 )Z 2 + T 2 S 2 Z 3 (7.15) 

Again we shall interpret this relation as an elliptic curve with Neron parameters S 
and T . The discriminant and j invariant for the curve are determined to be: 

A = -16(-4(S* 2 + T 2 ) 3 + 27T 4 S* 4 ) 

j = -16(-4 • 1728(5 2 + T 2 ) 3 /A ( ' 16) 

These examples suffice to illustrate the close relationship between the ground rings 
of (c = 1) string theory and algebraic curves. 

7.3 Modular forms and elliptic curves 

Consider the ordinary fermionic propagator in two Euclidean dimensions given 
by gi(z) = 2tt (0\ip(z)ip(0)\ 0) = - . This function satisfies the simple relation 
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{d 2 9l ) 2 = A{d gi f or 

(2ird z T zz ) 2 = 4(2nT zz ) 3 (7.17) 

with 2ttT zz = dg\ = On a two torus the relation is more complicated and is 
given by: 

(2ird z T zz f = 4(2ttT zz ) 3 - Q0G 4 (X)2itT zz - 120G 6 (A) (7.18) 
This equation is in one to one correspondence with the elliptic equation [43] 

Y 2 Z = 4X 3 - G 4 (A)XZ 2 - G 6 (A)Z 3 (7.19) 

Setting Z = 4 we obtain an elliptic equation with a = —4G 4 (A) and b = —16^6 (A). 
The correspondence is through (X, Y) <-> (T zz ,d z T zz ) . On the two torus the stress 
energy tensor can be expressed in terms of the Weierstrass elliptic functions through 



2^) = -, + E( I _ F -- 5 ) (7.20) 



1 , ! I 

uieA 

with the lattice functions 

G 4 (A) = 60 E A 
G 6 (A) = 140 £ 

Using the modular parameter of the torus the lattice vectors are of the form uo = 
2muj\ + 2nuj2 = 2ixR\m + 2tx{t\ + iR 2 )n = 27ri? 1 (m + nr) .and the above expres- 
sions define Eisenstein modular forms. The discriminant and j invariant also define 
modular forms and are given by: 



A = G 4 (A) 3 - 27G 6 (A) 2 = q f[ (1 - q n ) 24 = 2~ 8 6' 1 S = 2- 8 0f0f0f (7.22) 

(7.23) 



71=1 

1728(4a 3 ) 1728G 2 (A) 3 



J 4a 3 + 276 2 G 4 (A) 3 -27G 6 (A) 2 
The j function is a modular form of weight [43] and can be written as: 

. (9! + el + elf 



(7.24) 



8A 

This review is enough to illustrate the close connection between modular forms and 
elliptic equations. 

7.4 Dirichlet series, modular forms and elliptic equations 

The correspondence between theta and modular functions and Dirichlet series is 
transparent through the Mellin transform. Consider the formula for the Riemann 
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and Lerch zeta functions: 



S oo 

CM = 21-^-1 2^-1 ^"^ / v rS/2(04(o ^ r) + 02(o|ir) - 03(o|ir)) (7 ' 25) 

o 



s/2 °° 

<f>(x, s, a) + e" 2 ™ </>(-£, s, 1 - a) = ^— — / drr 5 / 2 ' 1 / 2 £ e -<n+ a fr ^inx 

neZ 

(7.26) 

A further example is the Ramanujan zeta function which in this case is the Mellin 
transform of the modular A function: 

oo 

F(s) = £ T(n)n- S 

71=1 

(7.27) 



r S 2 i v ; 



Where the r(n) coefficients are defined from the series expansion 

A = £ r(n)e _27r . We treat this and other zeta functions more fully in the next 

n=0 

section when we discuss string modifications to the fermionic Green's function. 

Another way of relating algebraic equations to modular functions is through 
elliptic genera. In this case the algebraic equation defines a manifold and the path 
integral of a sigma model with that manifold as a target space defines an invariant 
of that space. Computing the torus partition function then yields the relation to 
modular functions . An example of this type of analysis for singular K% and Calabi- 
Yau manifolds is given in [48]. 

Finally when Dirichlet series have Euler products representations and are asso- 
ciated with modular forms there is a direct relationship between elliptic curves E 
over the rationals and Dirichlet series. In this case the Euler product is directly 
associated with the elliptic curve by forming [43]: 

F(E,s)=T\z t— Th " (7-28) 

with 

a p =p+l-N p , 1,-1,0 (7.29) 

if p is good, split nodal, nonsplit nodal, cuspidal respectively and where N p is the 
number of points (x, y) satisfying the elliptic equation E modulo a prime p. An 
example is the Ramanujan zeta function defined above which has the Euler product 
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representation: 



1 3 

F(s) = J] (1 - r(p)p- s +p 11 - 2s )- 1 :(*>-) (7.30) 



This can be compared with the relatively simple Euler product for the Riemann 
zeta function 

^) = T\ T ^— s =ZfZ b (7.31) 
In the above we can express the product over primes as two factors 

71=2 

and 

oo 

Zf= II (1-M- S ) (7-33) 

m<r=2 

which have the interpretation of bonic and fermion partition functions with Hamil- 
tonian 

oo oo 

H = ( lo g( m ) + log(r))6_ m - r b m ,r + X! (log(n)a_„a„ (7.34) 

m<r=2 n=2 

And single particle energies 

e m>r = log(m) + log(r) 
e n = log(n) 

The reason for writing the zeta function this way is that if n = mr for some (m, r) 
than the factors in the numerator and denominator cancel and one is left with the 
product over primes. 

8 Other representations of the fermionic Green's func- 
tion 

The basic representation of the fermionic Green's function we have been using is the 
mode expansion representation given by: 

G(x, t; x', t') = (0\i>(x, t)^(x', 0| 0) = £ Y> n (x, (8-1) 

n 

where ip n (x,t) solves the generalized Dirac equation. However many other repre- 
sentations of the fermionic Green's function can be derived including the energy- 
momentum representation, point particle path integral representation [49], [50], [51], 
[52], random walk representation [53], [54], spin model representation [55], ID Ising 
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model representation [56], [57], [58], 2D Ising model representation [59] etc. As we 
have shown the zeta function can be interpreted as a Green's function so we should 
be able to obtain many of the representations for the zeta function using these 
methods. In this section we explore some of the ways of representing the fermionic 
Green's function associated with the Riemann, Lerch and Ramanujan zeta functions. 

8.1 (E,p) Energy momentum integral representation 
Usual case: linear dispersion E = p 

Recall that the usual massless 2d Dirac equation is given by 

(l°dt + r y 1 d x )ip(x, t) =0 (8.2) 

The Fourier transform of this equation is then 

( 1 °E + 1 1 p)4>(p,E)=0 (8.3) 

The Green's function G is then the inverse of the differential operator in (8.2) and 
is expressed as 

G(x, t; 0, 0) = JL J dpdE ^ WP^ -iEHtp* (fU) 
Writing out the gamma matrices the Green's function is a 2x2 matrix given by: 



f If e -iEt+ipx ^ 



G(M;0,0) 



1 f e~ 1 ' 

wi dpdE — 



(2tt) 2 J * p-E + ie 
dpdE 



.5) 



V (2vr) 2 J * p + E + ie J 



We shall concentrate on the upper right part of the matrix associated with the 
propagation of a right moving particle. Then the integral becomes 

G 12 (M;0,0) = — / dpdE L^ e -iEt + i P x = (8 6) 

(2tt) 2 J p-E + is 2-rrx-t 

in agreement with (3.5). 

Nonstandard case : Logarithmic dispersion aE = \og(ap + 1) 
Now applying this representation to the Dirac equation: 

(l°d t + 7 1 ^log(a^ x + l)1>(x,t) = (8.7) 
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The Fourier transform form of the equation is 

( 7 °£ + 7 1 - log(ap + l)^(p, E) = (8.8) 

Inverting this operator and again picking out the right moving portion the modified 
form of (8.8) becomes 

G 12 (x, t; 0, 0) = — ^ / dpdE, e -iEt+i P x (8 9) 

W J - log(ap + 1) - E + is 

If x is periodic with radius R the momentum is quantized so the p integral becomes 
a sum and we have: 

g 12 (m;0,o) = t^t^E / ^ e -^*+w* (8 . 10) 

a it 

Then if we assume R = a which is the case relevant to the Riemann zeta function 
we have: 

Gi 2 (M;0,0) = 7^T2- E / dE l " e - iEt+in */ a (8.11) 

a n=1 J log(n) _ ^ + i£ 

Performing the contour integral over E we obtain the ip function 

, OO -i 

G(x,t;0,0)| R „ = *JL, ; ) = £ -^e^/- (8.12) 

n=l 

For x = we obtain the Riemann zeta function: 



Gi 2 (x = 0,t;0,0)| R=a = ic(^) (8.13) 
in agreement with the mode expansion representation. 

Nonstandard case: exponential dispersion ap + 1 = e aE 

On the other hand if we start with the second form of the modified Dirac equation: 

{l °h {eta9t ~ 1} + ^dxM^t) = (8-14) 
The Fourier transform form of the equation is 

tf±(e aE -l) + *y 1 P ) 1 (,(p,E)=0 (8.15) 
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Then the modified form of (8.6) becomes 

■^-g 3 >(x,s) = G 12 (x,t + ia;0,0) 

= _i_ / dpdE = e -iE(t + ia) +ipX (8.16) 

^ J p -I( e ^-l) +ie 

Again If x is periodic with radius R the momentum is quantized so the p integral 
becomes a sum and we have: 

<?4' (x, s) = G41 (x, t + ia; 0, 0) 

1 1 ^ r .„ 1 



d£ e -iE(t+ia)+inx/R (8.17) 

(2-) 2 i?^i "I (e «*_ 1)+ie 

it a 

x 

Now we can use the formula X) = 7 — - — with u = — (e a — 1) 

n=-oo u — n sm(mi) a 

and W = e aE to write: 

54' (x, s) = G4' (x, t + ia; 0, 0) 
R 

-in — (W— 1) -R, TIr lX X —it /o io\ 

1 1 fdWire a y 1 i-(W-l)-( + 1 (8-18) 

= / ^ e a RW a 

sin(7r— (W - 1)) 
-R 

Switching variables tou = 2iir — W we write this as: 

a 



<7 4 / (x, s) = G41 (x, t + ia; 0, 0) 



R 




1 27TI — 

= -e a e " 


-ix/a/R 

a 


a 






1 2-wi — 

= -e a e" 


-ix/ar 

a 


a 



- s + l -«(-— 



' ',/ 7 e ,_ e «// 27ri? (8-19) 

L (27r)S -l^( S -l)/2 r(1 _ s)0( ;R j _ flj x _ X 

a ZttR 

For the case relevant to the Riemann zeta function R = a and x = ira this simplifies 
considerably to: 

9a'(x = ira, s) = G#{x = ira, t + ia; 0, 0)\R= a 

= i(2^) s - 1 e i7r ( s - 1 )/ 2 r(l - s)(2^ s - 1)C(1 - s) (b ' 2{)) 

Using the functional equation for the zeta function this can be written as: 

54' (x = ira, s) = G^{x = ira, t + ia; 0, 0)|# =a 

111 , (8.21) 

= -- — (2 1 " 5 - l)C(s) 

a 2 1- e l7rsV AV y 



78 



Thus the integral (E, p) form of the Green's function leads to the similar conclusion 
as the mode sum: The physical transition probability of the fermion leaving from 
(x = 0, t = 0) and arriving to (x = ira, t) vanishes at the Riemann zeros. 



8.2 Symmetry properties of the Green's function 

The usual Dirac fermion is two dimensions has a number of symmetries including 
translation, Lorentz rotation, scaling and especially conformal transformations. To 
determine the symmetries of the fermionic Green's function that leads to the Rie- 
mann zeta function we perform a Wick rotation to Euclidean space which changes 
the Green's function to: 

j°-(e iaE -1)+^ 
G af3 (x, t; 0, 0) = / dpdE—^ — = e -iEt+i P x ( 822 ) 

P 2 - A,(e iaE - l) 2 
cr 

Upon Wick rotation the nonstandard Dirac equation becomes the difference equa- 
tion: 

^ l -(e^-l) +1 l dMx,t) ( ^ 23) 
= ( 7 °(V(x,t + a) -^(x^ + ^d^x^t) = 
which is a field equation from the discrete field theory: 

1 = J dtdx(4>(x,t)(~f°(ij(x,t + o) - ip(x,t)) + 7 1 <9 X )V(M~) (8-24) 

The Green's function is modified under scale transformation as: 

G(Xx, Xi; 0, 0; XL, Xa) = X^G{x, t; 0, 0; L, a) 

7 oI (e ^_l )+7 i p 
= A" 1 E JdE a e -iEt+ iP x (8.25) 

= 27m p 2 + _( e ai£_ 1) 2 

p L a 

with A an integer. . The representation of discrete rotation and scale invariance can 
be found be studying the transformation properties of the denominator in (8.22) 
which results from applying the discrete Dirac operator twice. Setting K(p, E) = 

p 2 7:(e aiE — l) 2 = k 2 T7(e taw — l) 2 the generators of discrete rotation invariance 

a z or 

are found to be £q — lo = — (e iaw — l)dk — kd u and dilatation 

ai 

+ = dkk + 9^ — (e iau} — 1). In terms of position coordinates these are given by 
ai 

Iq-Iq = x-(e adt -l)d k -id x and £ +£ = xd x + t-(e adt -1). Further understanding 

a a 
of the symmetry properties can be found by using discrete conformal field theory as 

in [60] or by using methods of stress energy tensors for higher derivative theories. 
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We hope to return to these issues in future work. 



8.3 Green's function from periodization 

In the usual case (E = p) we can form the Green's function on a cylinder from the 

Green's function on a plane by periodization. That is beginning with 

i 

2nGi(x, t; 0, 0) = we form C?2(x,i;0, 0) through: 



2^G 2 (M;0,0) = ^G 1 (x-2irRn,t;0,0) = ]T _ - (8.26) 

n=— oo n=— oo 

oo 1 ire~ iwu — 2iir 

Now using the formula J2 = 1 — r = ^-^ we have: 

n=-oo u — n Sm{TTU) 1 — e /7rat 

2nG 2 (x,fA0) = ^ i _ e l_ t)/R (8.27) 

in agreement with (3.16). 

We can use a similar method to construct the nonstandard Green's function 
G4(x,t;0,0) . Beginning with the Green's function gy (x, s) = 2-JTGy(x, t— ia; 0, 0) = 

-(— ix/a) s ~ l e~' lx l a T{l — s) associated with the dispersion relation E = — log(ap + l) 
a a 
and with s = it + a we form: 

oo 

54 (x, s) = 2^(74 (x, t — ia; 0, 0) = J2 93' ( x ~ 2irRn, s) 

n=— oo 

-. oo 

= r(l-s)- ^ (-i(x - 27rn J R)/a) s - 1 e-^ M )/ a 

an= "°° x i? (8-28) 

= r ( i-.)I(V I (-fi-^r^ 

a a it 

n=— oo 

= -5(-)V(o £ 5.*.-) 
it a 27rit a 

in agreement with (3.52). 

8.4 Point particle path integral representation 

Another way of representing the Dirac Green's function is as a point particle path 
integral from point (0, 0) to (x, t) . This is given by: 

G(x,t;0,0) = J DeoDxoDxDtDpDEDXe* 1 (8.29) 

In the path integral (eo,Xo) represent a (0+1) dimensional gravitational and Rarita 
Schwinger field, (£,£) represent worldline fermions which eventually become the 
space-time gamma matrices, (x,t,p,E) are the position, time, momentum and en- 



80 



ergy associated with the point particle. 

For the usual case of the Dirac particle the (0+1) point particle action is: 

1 = JdT(p±-Ei-^(p 2 -E 2 )- itt - i|| - ipXoZ ~ iExot) (8.30) 

This point particle action has (super) reparametrization symmetry associated 
with the (0 + 1) dimensional supergravity. As discussed in several references [49], 
[50], [51], [52] the Green's function reduces to an integral over a single Teichmuller 
parameter C = J eo(r)dr and the zero modes of p and E . In [19] the modular group 

ri 

was described by by the Z2 transform eo(r) = — eo(— r) . Under this transformation 
the modular parameter changes sign and one restricts it to a fundamental region 
F = [0, 00] . Upon Wick rotation the Euclidean Green's function is given by: 

G(x, t; 0,0)= J dCdpdEe~ Cp2 - cE2 (£p + lE)e ipx ~ iEl (8.31) 

F 

After integrating over C and identifying the zero modes of (£,£) with the gamma 
matrices we obtain the usual form of the Dirac Green's function. 

In our case we can follow essentially the same procedure except our point particle 
action is more complicated, also the spatial direction is periodic or antiperiodic which 
causes the momentum to be quantized. The point particle action is: 

/ = Jdripx -Ei- |V - \{e aE - If) 

n 1 1 2 a (8.32) 

a a 

The path integral reduces to modular integral as before and we have: 
2ttG(x, t; 0, 0) 

1 f -Cp 2 +c-^(e iaE - l) 2 ~l ■ p ■ m (8-33) 

= W~ dCdpdEe a 2K ' + £-(e mE - l)yP x ~ lEt V ' 

27T J a 
F 

The zero modes of the fields become the gamma matrices 7 1 ,7° so after inte- 
grating over C we obtain the Green's function (8.22). To obtain the theta function 
representation of the Green's function we send p = — , R = ma and w = e iaE the 
Green's function reduces to: 

G(x, i; 0, 0) 

n 2 inx 1 2 

= ±\J dC^e'^ + TdWe- C ^ iW " 1} - +i(W - 1)W-^ 

p n 

(8.34) 
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For simplicity consider the case x = and R = a which is relevant to the Riemann 
zeta function. 



F W 

Now use the fact that 



g 4 > (0, s) = 2vrG 4 ' (0, t + ia; 0, 0) 

- fdC^e- Cn2 e- c W-V\n - 1 + W)W' S+1 ^ 



g fdCe- Cn2 e- c ( w -V 2 (n + W-l) 

30 

g JdCe- c( - n +V 2 e- cw2 (n + W + l) = £ / dCe~ Cn2 e~ cw V 



n=— oo 

oo 



(8.36) 



We have: 



54 



(0, s) = 2ttG 4 , (0, t + ia; 0, 0) = J dC^- ^ e" Cn2 e" CH/2 W~ s+2 (8.37) 



Performing the integral over W we have the theta representation. In terms of the 
zeta function this becomes: 

g 4 ,(0,s) = 27rG 4 '(<M + »a;0,0) = - — - COO (8-38) 

a sin(7rs/zj 

8.5 String modified Green's function 

The representation of the zeta function using theta functions represents a convenient 
way to study string modifications to the Green's function. This is because of the 
form of the expression : 

1 r - c {n + A? 
2ttG(x = 0, t; 0, 0) = — — / dC(^ e R 2 ) e -CM* c -u/a (8 39) 

where the modified dispersion relation including a mass M is 

( I e a E) 2 = (!L + i)2 + M 2 (8.40) 

a R 

or E 1 = — log(a(( — - — ) 2 + M 2 )). For the choice A = — this becomes: E = 
la R a 

— log((a— + l) 2 + a 2 M 2 )) which reduces to the usual dispersion relationE 1 = 
2a R 

I Tl 

— log(a— + 1) if M = 0. To study the string modification to the Green's func- 
a R 

tion from a point like configuration to another point like configuration one replaces 
the e~ CM2 factor with tr(e~ CM2 ) where a'M 2 = N and N is the occupation number 
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in the string oscillator description. Expanding the trace as 

oo 

tr(e- CM2 ) = ]T p(N)e- CN / a ' (8.41) 



N=l 

The Green's function becomes: 



2ttG(x = 0,t;0,0) 

-, „ oo X n + A ) 2 oo 

= fW2)/^S e R2 )E P We- CN/a 'C-* t/2a (8 .42) 



F 



n=l,N=l ti a 

Two dimensional heterotic string 

The actual value of the degeneracy factor p(N) depends on the string theory con- 
sidered. For example the choice: 

oo /) 1 2 oo oo 

tr(e- CM2 ) =,11(1 + <f ) 24 = 2" 8 ^ = £ p 12 (7V)^ = £ r'(A^ (8.43) 

n=l 1 N=l N=l 

where q = e~ c ^ a describes a sector in the fermionic formulation of Eg x SO (8) or 
50(24) heterotic two dimensional string theory [61] and the subscript denotes the 
dimension of the root lattice of the group. Taking the leading term in the series 
over n by sending R — > and setting a 2 = a' we have: 



oo 



G(x = 0, t; 0, 0)| R=0 = ]T Pi2(N)N- u / 2a = ]T r f (N)N- u / 2a (8.44) 

JV=1 N=l 

Ramanujan Zeta function 

Another choice is 

oo oo oo 

tr (e- CM *) = q f[ (1 - q n ) 24 = 2" 8 ^ 8 = £ p(N)q N = ]T T(N)q N (8.45) 

n=l N=l N=l 

This form is related to the Ramanujan zeta function [7]. 

oo 

F(s) = T{n)n~ s (8.46) 

n=l 

The Ramanujan zeta function has similar properties to the Riemann zeta function 
with functional relation: 

(2tt)- s T(s)F(s) = (2vr)-( 12 - s )r(12 - s)F(12 - s) (8.47) 
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The analogue of the Riemann hypothesis for the Ramanujan zeta function is the 

12 13 

absence of zeros in the region — < Res < — . We can use the formula 

oo 

Tr(e- CM2 ) = q ]J (1 - q n ) 2A = 2- 8 6[ 8 (8.48) 
m=l 

where q = e 2l7TT to write the Green's function in theta form : 
g 5 (x, s) = 2irG 5 (x, t; 0, ia) 




(8.49) 



We can also introduce the string generalization of the states 



oo 

v')= J2 e- aE "> N / 2 (-l) n \n,N) 




n,N) 



(8.50) 



Then the amplitude for transition from the state v to v' is given by: 



oo 




(8.51) 



n,N=l 



Using the dispersion relation (8.40) we have: 




(8.52) 



The transition amplitude is: 




7l,JV=l 



(8.53) 



To compare with the Kaluza-Klein neutrino oscillation we note that 




(8.54) 
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Whereas the Kaluza Klein expansion is given by [12]: 



oo 

'I V I 

\v) KK = J2 i = K N ) ( 8 - 55 ) 

n=0,N An + Rnv) . 2n 

So that the form is the same when m = —, [iy = A/R, fiy = M and a = a. Again 
the time variation of this state is very different in our case because of the nontrivial 
logarithmic dispersion relation. 

Taking the limit of R — > only the first term contributes as the other terms are 
infinitely massive and we have the transition probability: 

00 a 2 

p(i/-n/,t) = l(i/|i/,*)l 2 = l (^w) (_i *~ ff)/2 MW)l 2 ( 8 - 56 ) 

n,N=l a 

Ten dimensional open fermionic string 

Finally the open ten dimensional fermionic string corresponds to the choice [62]: 



Tr(e-^ 2 ) = ^ = ^ = 4= E Po{N)e-™/<*' (8.57) 

n (1 - g n ) 8 n (i - i n/2 ) 4 1 

n=l n=l 

In this case the dependent Green's function is given in the (E,p) representation by: 
G(x,t;0,0) = 

JdEd 9 pdC( 7 °-(e aE - l)+ Pa i )e- iEt+ipx e~ C{ ~a ieaE ~ 1))2 + ^Tr^ ^) 

(8.58) 

As we are interested in short distance behavior of the Green's function we examine 
the region of the integrand where e aE » 1 . Also we assume x = so that: 



G(0,t;0,0) 



JdEd^-e aE e^ Et [ dCe^a {e ^ Y+p2 + i£) £ p (N) e - CN / a ' 

Qj J at r\ 



JdE 7 °-e aE e- iEt { dCe° { £ p {N)e~ CN / a ' ' C^ 2 

a J N=Q 



N=0 



(8.59) 
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To extract the N dependence redefine 2aE = 2aE — log(iV) and C = CN 



G(0,t;0,0) 

~ JdE 7 °-e aB e- iEt I dCe~ G{ ~~a {eaE)2 + i£) e -C/«'c-»/ 2 f] Po (N)N~to +4 

a J N=0 

(8.60) 

So that the short distance Green's function is again described by a Dirchlet ex- 
pansion. The quantization in this case is coming from the string mass condition 
M% = TV/a'. 



9 Conclusions 

In this paper we represented the Lerch zeta function as a Green's function associated 
with a two dimensional Dirac equation and derived the Dirichlet functions including 
the Riemann zeta function as special cases of the relation. We then formulated the 

Riemann hypothesis in physical terms associated with the transition amplitude and 

1 2isn 
dispersion relation E = — log(l+ap) . The momentum p was quantized — — because 

a Lj 
the spatial direction was a chosen to be a circle of circumference L . The circle in 

turn obeyed the special relation L = 2ira .Within this framework we obtained a 

particular representation of the Riemann hypothesis as a nonzero condition on the 

transition probability between two states that is analogous to the phenomena of 

neutrino mixing. 

We discussed the parameter a from several viewpoints including mixing param- 
eters, duality, susy breaking, zeros in path integrals and as an extra dimension. We 
also studied the thermal partition function and thermal Green's function associated 

with the dispersion relation and derived the thermal divergence at temperature — 

a 

. We explained the relation of the thermal partition function of this theory to the 
thermal partition function of Julia, Spector and Bowick [27], [28], [29]. We studied 
the fermionic string theory associated with a logarithmic oscillator, its matrix model 
and the relation of the periodization of the inverted oscillator to the Riemann zeros. 
We discussed the relationship between the ground rings of (c = 1) string theory, 
elliptic curves, modular forms and Dirichlet series. Finally we derived several rep- 
resentations of the fermionic Green's function including the (E,p) representation, 
point particle path integral representation, and string modifications. For the string 
modified Green's functions we derived Dirichlet series for the 2D heterotic string, 
the Ramanujan r function and the open ten dimensional fermionic string. Other 
physical pictures and methods can be obtained by permuting the physical variables 
(x,t,p,E) or searching for Yang-Lee zeros of a one particle thermal partition func- 
tion. One can also consider the Green's function of a nonstandard scalar particle. 
We discuss some of these in Appendix A , B and C. 



86 



A Permutations of x,t,p,E and different physical real- 
izations of the Green's function 



One can obtain different physical realizations of the two point function simply by 
permuting the dynamical variables (x,t,p,E) . 

In the interpretation we have been using in this paper the point particle path 
integral representation of the Green's function is: 



x(l)=x,t(l)=t 

G(x,t;0,0)= J DxDtDEDpDeoDxoDZDte-'J ' Ldr (A.l) 

i(0)=0,t(0)=0 



where 



L=px-Ei + ^-e (p 2 - -^(e aE - If) + ixoPt + iXo~(e aE ~ ~ ~ %Z (A.2) 



The classical equations associated with this choice are given by: 

Energy : E = — log(ap + 1) 
a 1 

momentum : p = —{e aE — 1) 
dE a c 

velocity : v = — = — - — = ce~ aE (A.3) 
op ap + 1 

position : x = 1 = e ct 

ap + 1 
ap + 1 „t?x 
time : t = — x = e aE - 

c c 

Finally we assumed periodic boundary conditions in the x direction so that momen- 
i n 

turn was quantized p = — . 

R 

Sending x <-> t;p <-> E 

This is a simple permutation exchanging space with time and momentum with en- 

l> t x R 

ergy. The Green's function in this case becomes G(t, x; 0, 0) = ( — Y x ^ a <j)( , — , — ) 

a 2-itR a a 

The bosonic piece of the point particle Lagrangian is: 



L = Et — px — ^(\{e ap - If - E 2 ) (A.4) 
— ■ a 
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The classical equations associated with the Lagrangian are: 



p = - log(aE + 1) 
E= -{e ap - 1) 

t = — x = e~ ap cx 

aE + l 

aE + 1 t 

x = 1 = e ap - 

c c 

Tl 

In this case one takes time periodic so that E = — . 

R 

Sending x <-> p;t <-> E 

R E v R 
The Green's function in this case is G(p,E;0,0) = (—) ip/a 4>( , -, — ). 

a 2ttR a a 

The bosonic piece of the point particle Lagrangian is: 



L = -px + Ei- — (4r(e at - l) 2 -x 2 ) (A.6) 
2 a z 

The classical equations become 



ce at (A. 7) 



The trajectory 



t = - log(ax + 1) 

x=\e at -1) 

_ %E _ ax + 1 

dp c 

p = — - — E = e~ at cE 
ax + 1 

E = ^±lp = e^ 

c c 



x = \ e at - 1) (A.8) 
a 



comes from the HamiltonianS = axp + p which reduces to the Berry-Keating oper- 
ator E = apx at large x [2]. The trajectory is also the geodesic associated with the 

metric ds 2 = —dt 2 + e~ 2at dx 2 . In this case space is quantized as x = — and a and 

R 

R have units of energy. 



Sending x <-> E;t <-> p 

Finally sending x ^ E;t ^ p the Green's function becomes 
R E v R 

G(E,p;0,0) = (—Y E / a (j)( — -,-,—). The bosonic piece of the point particle La- 
a 2irR a a 



grangian is: 



L = -Ei + px-—(\{e ax -l) 2 -t 2 ) (A.9) 
2 a 1 
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The classical equations become: 

x = - log(ai + 1) 
f 

t = -(e ax - 1) 

v = ^- = = ce~ ax (A. 10) 

dp at + I K J 

E = —^—p = e~ ax cp 
at + 1 

at + 1 a:E -E 
c c 

In this case the trajectory is the null geodesic associated with the metric ds 2 = 

—e~ 2ax dt 2 + dx 2 . Here time is quantized as t = — and a and R have units of 

ii 

energy. 



B One particle partition function 

The one particle partition function associated with the dispersion relation E 
1 27rn 

— log(a— h 1) is given by: 

a Li 



m = t .-«. = f ■ , - + (-) (B ,, 



L 

This reduces to the Riemann zeta function for the special case = 1. 

2ira 

Because of the form of the dispersion relation the one particle partition function 
is equivalent to the partition function of a nonstandard quantum oscillator with 

energy frequency relation E = -log(a(n + -)u) + 1). Note this relation reduces to 

_^ a £ 

the usual relation E = (n+ -)u in the limit as a goes to zero. The partition function 
of the nonstandard oscillator is: 



4(3, li) = E e~^-^ n 

n=0 

= £ (a(n + i)w + l)-^ e -^ = (au,)-^/^, ^ 1 + — ) 
n= o ^ a z au 



(B.2) 



where we have introduced a chemical potential associated with the oscillator number 
operator. Continuing this function to the complex plane we see that when the special 
relation auj = 2 is obeyed as well as (3/j, = iri the Yang-Lee zeros of the partition 
function coincide with the zeta function zeros. The Green's function nonstandard 
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oscillator can be denned by the path integral: 

x(t)=x 

G(x,t;0,0) = J DxDpe^ dTL (B.3) 

x(0)=0 

1 1 p 2 

with L = px log(a-( h muj 2 x 2 ) + 1). The classical equations associated with 

a 2 m 
the nonstandard logarithmic oscillator are: 

E = - log(-(— + muj 2 x 2 ) + 1) 
a 2 m 

p/m p 

x = 



|(^w* 2 ) + i mea£ 

p = - = -gme^V** = 

- ( h x ) + 1 

2 2 m 
P I 2 2 ^ / a£ r \ 

h mw ? = -(e — 1) 

m a 

x = x cos(oje" aE t) + — sm(uje~ aE t) 
p = Po cos(uje~ aE t) — XQU)msm{uje~ aE t) 

Again in the nonstandard case the period depends on the energy and hence the 

2ir 2tt 2tt 1 

amplitude through T = = = — e aE = — (1 H — mu> 2 x\a) for the case of 

u)e~ atj uj uj 2 

po = 0. The quantization yields the same eigenfunctions as the usual oscillator 

because of the expansion: 

H^ n = E n ^ n = {-a) k -\u{h + -) k )^ n (B.5) 
k=i 1 

Finally although the partition function of the nonstandard oscillator is the same 

as the one particle Weyl fermion partition function with dispersion relation E = 
1 2im 

— log(a— h 1) the underlying theories are different. To see this note that the 

(i Li 

Green's function for the nonstandard oscillator is: 

OO OO -. 

G(x, t; 0, 0) = e- iEnt H n (x)H n (0) = £ (a(n + -)uj + iy it l a H n {x)H n {0) 

n=0 n=0 

(B.6) 

Whereas that for the nonstandard Weyl fermion the Green's function is: 

OO OO r) 

G(x, t; 0,0) = J2 e- lEnt e 2 * in ' L = £ (a=^ + i)~^ e ^/L (R7) 

n=0 n=0 

These are not the same although the partition function defined from the Green's 

2ir 

function z((3) = J dxG{x, i/3; x, 0) are identical if one maps uj to — . Note the con- 
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dition to obtain the Riemann zeta function aw = 2 is far from the limit to obtain 
the standard oscillator au <C 1. 



C Scalar field 

Green's function 

Most of the discussion has been with a fermion field obeying the dispersion relation 

p = —(e aE — 1). Many of the considerations apply to a scalar field &(x,t) as well. 
a 

Also the theory of a single Weyl fermion has a gravitational anomaly. Adding left 
moving scalar matter to the theory can cancel this anomaly as happens in the Type 
OB string model. Taking the square of the Dirac operator we obtain the nonstandard 
Klein Gordon equation 

(d 2 x + ^(e mdt -l) 2 )$(x,t) =0 (C.l) 

Tl 

On a circle momentum is quantized as p = — . Forming a series expansion for the 
scalar field: 

*) = E -m^onM*, t) (C2) 

Where a_ n and a n are creation and annihilation operators and ip n are the mode 

n 

i — X 

solutions ip n (x,t) = e n e R . The one particle energies are: 

E n = -\og(a- + l) (C.3) 
a R 

The scalar field expansion is given by: 



x 

*(*, t) = £ . 1 Un + ir^e n R + c.c (C.4) 



«^ 2 _ log(a _ +1) 

The bosonic or scalar Green's function is denoted by K and is expressed as: 

n 

oo ( a - + l)- lt / a e R 
K(x, t; 0, 0) = <0|$(x, t)$(0, 0)| 0) = 2 „ (C5) 

n=o 2-log(a— + 1) 
a R 
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The two point function between the scalar field and its time derivative is somewhat 
simpler and can be expressed in terms of the Lerch zeta function: 



(OR*™, 0)| 0) = i E (o^ + D-V-e*^ = 1 (f )^-^. 1, |) 

n=0 

(C.6) 



Statistical Mechanics 

The statistical mechanics of the scalar field is determined by the Bose-Einstein 
distribution with the multiparticle Hamiltonian: 

oo 

H = E £ n a-na n (C.7) 

n=0 

and the one particle energies ae n = log(l + an/R). The partition function is given 
by: 

oo oo 

Z = tr(e~ pH ) = JJ (1 - e-fc")- 1 = JJ (1 - (a - + l)~^/ a ) _1 (C.8) 

n=0 n=0 

As in the Fermi-Dirac case one can develop a series representation for the free energy. 
Taking the logarithm of (C.8) we have: 

oo 

log Z = - ]T log(l - (o^ + (C.9) 

n=0 

Now we can expand the logarithm appearing in each term of the sum using the 

OO y m 

formula — log(l — y) = — to obtain the series representation: 

i°g^ = E E 1 (4 + ir /3m/a = E 1 (-) /3m/a CH(^,-) (c.io) 

n=0m=l m=l 

For the special case R = a relevant to the Riemann zeta function we have: 

oo n 

log Z(R = a) = E -C(— ) (C.11) 

m=l 

The thermal divergence occurs at the critical inverse temperatures [3 C = — with m 

m 

an integer. 

The multiparticle density of states can be obtained from the inverse Laplace 
transform of the partition function as 

ioo+c 

a{E) = J dpZ(P)e^ E (C.12) 

— ioo+c 
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For the special case R = a the total energy of the multiparticle bosonic system is 
given by: 



Etot = ~J2 N n l °S 



n 



(C.13) 



71=2 



Where N n is the number of particles of momentum n . As in the Fermi-Dirac case 
the quantity cr(E) has interpretation in terms of multiplicative number theory. <r(E) 
is the number of ways of writing e aE as a product irrespective of order and allowing 
repeat factors. 

For example we can enumerate the states that yield total energy E = - log(120) 

a 



120 = 5! = 

5 • 3 • 2 • 2 • 2;a_ 5 a_ 3 a 3 _ 2 | 0) 
15 • 2 • 2 • 2;a_i 5 ai 2 | 0) 

10-3-2-2;a_i a- 3 a^ 2 l 0) 

6 • 5 • 2 • 2; a_ 6 a_ 5 a 2 2 \ 0) 

5 • 4 • 3 • 2; a_5a_4a_3a_2| 0) 



30 



20 
15 
12 
10 



2;a- 30 a 2 _ 2 \ 0) 



2;a_2oa-3a-2| 0} 
2; 0-150,-40-2! 0) 
2; a_i 2 a-5a-2| 0) 
3; a_ioa_ 4 a_3| 0) 
8-5-3; a_8O.-5a.-3l 0) 
6-5-4; a_6a_5a_4| 0) 



2; a_6oO-2 
3; a_4oa_3 
4;a_3 a_4 
5;a_240-5 
6;a_2oa-6 
8; a_i50_8 



60 
40 
30 
24 
20 
15 
12 

120; 0-120 1 0) 



0) 
0) 
0) 
0) 
0) 
0) 



10;a-i 2 a_i | 0) 



(C.14) 



Thus a (— — — -) = 20 separate states. The bosonic scalar field has more states than 
the fermionic case as it allows repeats. 
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